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In  this  report,  research  topics  in  area  of^spread  spectrum  communications  that  are, 

fundamental  to  the  successful  operation  of  an  overall  system  are  discussed.  I Z  wi 3 1  A+. — - 

- — seen— the-t  system  performance  can  be  significantly  enhanced  by  the  incorporation  of  the 

transform  domain  processing  techniques  described  below.  The  following  key  topics  have  been 
studied.  .  f 

i)  We__iiaALe — impie-mentreH  a  transform  domain  processing  intercept  receiver,  for  detection 
of  direct  sequence  BPSK  signals  in  the  presence  of  finite  bandwidth  interference  using  an 
adaptive  interference  rejection  technique.  The  receiver  performs  a  real-time  Fourier 
transformation  followed  by  one  of  two  interference  rejection  techniquesi^  (1)  Interference 
detection  and  excision  by  notch  filtering,  or  (2)  Interference  attenuation  using  soft- 
limiting.  Interference  rejection  is  followed  by  inverse  transformation  and  then  detection 
by  a  total  power  radiometer.  While  the  exper lineni-di  Lesuits  show  that  both  techniques  can 
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yield  significant  improvement  in  system  performance,  the  soft-limiting  scheme  appears  to 
outperform  the  notch  filter  technique. 


'T-  ii)  During  the  process  of  acquisition,  a  spread  spectrum  system  is  at  its  most 
vulnerable  stage,  since  it  has  no  processing  gain  working  in  its  favor.  We  have 
investigated  the  use  of  a  narrowband  interference  suppression  f iltejt^based  upon  transform 
domain  processing,  to  enhance  the  perf ormance— of— a- serial- acquisition  scheme  for  a 
direct  sequence  spread  spectrum  receiver L?  An  analytical  expression  for  the  probability 
of  error  has  been  obtained  for  the  case  when  a  sinusoidal  interferer  is  present,  and 
numerical  results  have  been  generated  to  illustrate  the  sensitivity  of  the  system  to 
various  parameters.  )  We  have  shown  that  the  probability  of  the  system  locking  on  to  an 
incorrect  pha^a— position  is  significantly  lower  when  the  suppression  filter  is  used  than 
it  is-ih  the  absence  of  the  filter. 


'w>  iii)  Some  experimental  work  has  been  performed  in  conjunction  with  a  novel  synchron¬ 
ization  and  PN  code  acquisition  technique  employing  transform  domain  processing  for  a 
DS  spread  spectrum  communication  receiver.  The  most  important  advantage  of  the  proposed 
scheme  is  the  insensitivity  to  the  delay  between  the  reference  PN  code  and  the  received 
signal,  in  addition  to  the  capability  of  adaptive  narrow-band  interference  rejection  and 
carrier  synchronization  using  the  cross-correlation  between  the  received  signal  and  the 
reference  PN  code.  -  In  addition,  data  demodulation  can  also  be  performed.  Thus,  it  is 
envisioned  that  j^eompletely  functional  spread  spectrum  receiver  can  be  implemented  in 
the  transf_orm~dbmain  v/ith  superior  narrowband  excision  capability. 

iv)  '  . ‘ 
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"We  have  also  initiated  research  for  applying  spread  spectrum  techniques  to 
optical  communication.  Transform  domain  processing  in  the  optical  domain  was  investi¬ 
gated  and  a  design  for  spectral  coding  at  optical  frequencies  has  been  developed.  The 
encoding  system  is  based  on  well  established  principles  of  optical  pulse  compression 
and  shaping.  Several  possible  implementations  of  this  system  have  been  considered,  along 
with  its  applications  to  fiber  optics,  laser  radar,  free  space  optical  communications  and 
other  systems.  ,  ' 
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APPLICATIONS  OF  TRANSFORM  DOMAIN  PROCESSING 


TO  THE  RECEPTION  AND  DETECTION  OF  SPREAD 
SPECTRUM  SIGNALS 


1.  INTRODUCTION 

The  intent  of  this  report  is  to  describe  research  in  some  of  the  more  important  areas 
of  spread  spectrum  communication.  The  specific  topics  include  intercept  receiver,  acquisi¬ 
tion  of  direct  sequence  spread  spectrum  signal  and  spread  spectrum  techniques  in  optica] 
communication.  The  concept  that  ties  these  seemingly  disjoint  areas  together  is  the  use 
of  so-called  transform  domain  processing  which  is  discussed  in  the  next  section.  This  is 
followed  by  the  report  on  the  actual  research  performed.  Appendix  A  lists  the  papers  pub¬ 
lished  and  to  be  published  under  the  auspices  of  this  contract.  A  list  of  papers  presented 
in  national  and  international  meetings  is  included  as  Appendix  B. 

2.  TRANSFORM  DOMAIN  PROCESSING 

The  technique  of  transform  domain  processing  is  described  in  [l] .  In  this  technique,  a 
tapped  delay  line,  typically  implemented  with  a  surface  acoustic  wave  (SAW)  device  that 
has  a  chirp  impulse  response  built  into  the  taps,  is  used  as  a  real-time  Fourier  transformer. 
As  described  fully  in  [1],  a  notch-  filter  is  easily  implemented  by  Fourier  transforming  the 
received  waveform,  using  an  on-off  switch  to  perform  the  notching  operation,  and  then 
inverse  transforming.  A  block  diagram  of  this  system  is  shown  in  Figure  1. 

It  is  clear  that  the  heart  of  this  scheme  is  the  ability  to  perform  a  real-  time  Fourier 
transform,  and  the  technique  to  do  so  has  been  described  in  many  references,  (e.g.  [1] 
and  [3]).  Briefly,  the  Fourier  transformation  is  accomplished  in  real-time  by  a  SAW  device 
in  the  following  manner.  If  a  signal  f(t)cos(cu0£  +  0t2)  (that  is,  a  signal  f(t)  modulating 
a  linear  FM  or  chirp  waveform)  is  convolved  with  the  signal  cos (u>at  —  pt2),  the  result  of 
that  convolution  will  be  the  Fourier  transform  of  f(t).  That  is,  the  output  of,  say,  a  SAW 
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convolver  with  the  above  two  signals  as  its  two  inputs  (recall  a  SAW  convolver  is  an  active 
device  requiring  two  inputs)  will  be  a  bandpass  waveform  whose  inphase  and  quadrature 
components  correspond  to  the  real  and  imaginary  parts  of  F(w),  the  transform  of  f{t), 
respectively.  In  particular,  the  output  will  be  F(cu)  over  the  range  Loe[2{$T,2f}A),  where 
T  is  the  symbol  duration  and  A  is  the  time  for  the  acoustic  wave  to  travel  the  length  of 
the  device.  Alternatively,  rather  than  use  the  convolver,  if  one  implements  a  tapped  delay 
line  with  tap  coefficients  given  by  samples  of  the  unmodulated  chirp  signal  spaced  (n/uc) 
seconds  apart,  where  u>c  is  the  bandwidth  of  the  chirp  waveform,  one  can  obtain  F(w)  as 
the  output  of  the  delay  line  when  f(t)cos(u;0<  +  /3t 2)  is  the  input,  in  either  case,  the  radian 
frequency  variable  u>  will  be  a  linear  function  of  time,  so  that  F(u>)  will  be  generated  at 
the  device  output  in  real-time. 

A  detailed  performance  analysis  of  a  transform  domain  processing  system  is  presented 
in  [l]  and  will  not  be  duplicated  here.  Rather,  recently  obtained  experimental  results  ([4]) 
illustrating  the  behavior  of  the  system  will  be  briefly  discussed  in  order  to  provide  a 
perspective  on  how  well  the  system  can  perform. 

The  system  actually  implemented  is  shown  in  Figure  2,  and  it  can  be  shown  in  a 
straightforward  manner  that  the  configuration  is  equivalent  to  that  of  Figure  1.  The 
spread  spectrum  code  used  for  the  experiment  is  a  63  chip  PN  sequence  with  a  chip  rate 
of  1.875  MHz  and  is  shown  in  Figure  3(a).  Figure  3(b)  shows  the  Fourier  transform  of 
this  signal  as  obtained  at  the  output  of  a  real-time  Fourier  transformer  which  uses  chirp 
devices  with  center  frequencies  of  15  MHz,  bandwidth  of  7  MHz,  interaction  times  of  117 
/isec  and  chirp  rates  of  (3  =  3  x  1010  Hz/sec.  Thus,  in  the  cu-domain,  1  psec  corresponds 
to  60  KHz  and  the  main  lobe  of  S (cu)  has  a  width  of  3.75  MHz.  Figure  3(c)  shows  the 
Fourier  transform  of  a  single  tone  interferer  gated  by  a  pulse  of  duration  33.6  p sec.  The 
width  of  the  center  lobe  corresponds  to  60  KHz. 

Some  qualitative  results  of  processing  the  received  waveform  are  shown  in  Figure  4. 


3 


Figure  4(a)  shows  the  output  of  the  correlator  in  the  absence  of  any  noise,  and  Figure 
4(b)  shows  the  unfiltered  correlator  output  pattern  corresponding  to  an  interference  power 
which  is  20  db  larger  than  the  signal  power.  Figure  4(c)  shows  the  same  pattern  when  the 
transform  domain  filter  with  a  notch  width  corresponding  to  180  KHz  is  switched  on.  The 
notched  spectrum,  R(w)H(w),  is  shown  in  Figure  4(d). 

Figure  5  presents  curves  of  probability  of  error  vs  energy  per  bit  to  noise  spectral 
density  ratio  (E b/v)  for  the  case  of  single  tone  interference.  The  curves  are  parametrically 
normalized  by  the  offset  frequency,  6uo,  from  the  center  and  by  whether  or  not  the  notch 
is  employed.  When  the  notch  is  indeed  used,  the  notch  width,  Su,  is  fixed  at  67r.  For  Sui o, 
7T  corresponds  to  30  KHz.  Theoretical  curves  are  also  presented  in  the  same  figure  and  the 
agreement  is  within  a  fraction  of  a  dB.  To  achieve  adequate  phase  averaging,  the  single 
tone  interferer  is  phase  modulated  with  a  phase  excursion  of  ±n  radian  at  100  Hz.  The 
signal  power  to  interference  power  ratio  is  -20  dB  for  all  the  measurements  reported  here. 

Figure  6  shows  curves  of  probability  of  error  vs  notch  width  for  two  offset  frequencies, 
<5u;o  =  0.57T  and  Su 0  =  327r,  for  E^/p  =  10  dB.  It  is  seen  that  the  measured  curves  for 
Pe  fall  between  the  theoretical  curves  for  E#/t?  =  9  dB  and  10  dB.  For  the  case  of  <5w  = 
0.57T,  an  optimum  notch  width  around  107r  has  been  found.  However,  for  Suo  =  327 r,  the 
interference  sidelobes  do  not  overlap  much  with  signal  and  the  Pe  curve  is  more  or  less 
independent  of  the  notch  width  beyond  a  certain  value. 

Figure  7  presents  results  corresponding  to  the  interferer  being  phase  modulated  with 
a  phase  excursion  of  and  a  phase  modulating  frequency  of  im  =  0.5  KHz.  Comparing 
curve  (a)  of  Figure  7  with  curve  (d)  in  Figure  5,  it  is  seen  that  no  appreciable  difference 
is  found.  Thus  the  adaptive  system  of  tracking  the  center  frequency  of  the  interferer  is 
effective  for  this  case.  If  we  increase  to  5  KHz,  curve  (b)  is  obtained,  and  it  is  now  seen 
that  significant  degradation  has  occurred.  However,  if  we  double  the  notch  width  to  127t, 
curve  (c)  is  obtained,  which  again  shows  the  interference  suppression  is  effective. 
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In  an  actual  system,  the  center  frequency  is  unknown  and  probably  varying  with 
time.  Thus  a  more  meaningful  measurement  corresponds  to  allowing  the  center  frequency 
to  be  swept  slowly.  This  is  shown  in  Figure  8  for  a  swept  frequency  of  1  sec/sweep  and 
a  frequency  excursion  of  ±1  MHz.  It  is  seen  that  the  interference  suppression  is  quite 
effective  and  is  not  a  very  sensitive  function  of  notch  width  over  a  broad  range  of  values. 

Finally,  we  have  also  studied  the  effect  of  amplitude  weighting  the  received  signal  in  a 
spread  spectrum  transform  domain  processing  receiver.  Both  experimental  and  theoretical 
results  of  probability  of  error  have  been  generated  for  direct  sequence  systems  using  a 
63  chip  spreading  sequence  with  non-contiguous  data.  The  DS  signal  was  received  in 
the  presence  of  narrowband  interference  and  additive  white  Gaussian  noise,  and  good 
agreement  between  theory  and  experiment  were  demonstrated.  Results  were  given  in  [8] 
for  both  weighted  and  unweighted  input  signals  at  various  notch  bandwidths. 

3-  ADAPTIVE  NARROW-BAND  INTERFERENCE  REJECTION 
IN  A  DS  SPREAD  SPECTRUM  INTERCEPT  RECEIVER  USING 
TRANSFORM  DOMAIN  SIGNAL  PROCESSING  TECHNIQUES 
The  experimental  results  are  discussed  in  Papers  #1,  #2  and  #9  in  Appendix  A. 
The  theoretical  results  are  discussed  in  paper  #10.  As  papers  #9  and  #10  are  yet  to  be 
published,  they  are  included  as  Appendix  C  and  D  respectively. 

4.  INTERFERENCE  SUPPRESSION  TO  AID  ACQUISITION  IN 

DIRECT  SEQUENCE  SPREAD  SPECTRUM  COMMUNICATIONS 

The  results  are  discussed  in  Papers  #3,  #4,  #5,  #6  and  #8  in  Appendix  A.  As  the 
papers  have  already  been  published,  we  only  include  the  following  abstracts. 

ABSTRACT  OF  #3 

During  the  process  of  acquisition,  a  spread  spectrum  system  is  at  its  most  vulnerable 
stage,  since  it  has  no  processing  gain  working  in  its  favor.  In  this  paper,  a  narrowband 
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interference  suppression  filter  is  used  to  enhance  the  performance  of  a  serial  search  ac¬ 
quisition  scheme  for  a  direct  sequence  spread  spectrum  receiver.  An  analytical  expression 
for  the  probability  of  error  is  presented,  and  numerical  results  are  used  to  illustrate  the 
sensitivity  to  various  system  parameters. 

ABSTRACT  OF  #4 

Spread  spectrum  communication  systems  have  many  applications,  including  inter¬ 
ference  rejection,  multiple  accessing,  multipath  suppression,  low  probability  of  intercept 
transmission  and  accurate  ranging.  Of  all  the  potential  applications,  the  ability  of  a  spread 
spectrum  system  to  withstand  interference,  both  intentional  and  unintentional,  is  proba¬ 
bly  its  greatest  asset.  Since  any  spectrum  receiver  can  only  suppress  a  given  amount  of 
interference,  if  the  level  of  interference  becomes  too  great,  the  system  will  not  function 
properly. 

Even  under  these  latter  circumstances,  however,  other  techniques,  which  enhance  the 
performance  of  the  system  over  and  above  the  performance  improvement  that  comes  au¬ 
tomatically  to  systems  simply  by  employing  spread  spectrum,  are  available  for  use.  These 
techniques  typically  involve  some  type  of  additional  signal  processing  and  are  the  subject 
of  this  paper.  In  particular,  two  general  types  of  narrow-band  interference  suppression 
schemes  are  discussed  in  depth,  those  based  upon  least-mean  square  estimation  techniques, 
and  those  based  upon  transform  domain  processing  structures. 

ABSTRACT  OF  #5 

Spread  spectrum  communication  systems  have  many  applications,  including  inter¬ 
ference  rejection,  multiple  accessing,  multi-path  suppression,  low  probability  of  intercept 
transmission,  and  accurate  ranging.  Of  all  the  potential  applications,  the  ability  of  a 
spread  spectrum  system  to  withstand  interference,  both  intentional  and  unintentional,  is 
probably  its  greatest  asset.  Of  course,  any  spread  spectrum  receiver  can  only  suppress  a 
given  amount  of  interference;  if  the  level  of  interference  becomes  too  great,  the  system  will 
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not  function  properly. 

Even  under  these  latter  circumstances,  however,  other  techniques,  which  enhance  the 
performance  of  the  system  over  and  above  the  performance  improvement  that  comes  au¬ 
tomatically  to  systems  simply  by  employing  spread  spectrum,  are  available  for  use.  These 
techniques  typically  involve  some  type  of  additional  signal  processing  and  are  the  subject 
of  this  paper.  In  particular,  two  general  types  of  narrow-band  interference  suppression 
schemes  are  discussed  in  depth,  and  a  short  overview  is  presented  for  several  other  tech¬ 
niques  as  w'ell.  The  two  classes  of  rejection  schemes  emphasized  in  the  paper  are  1)  those 
based  upon  least-mean  square  estimation  techniques,  and  2)  those  based  upon  transform 
domain  processing  structures. 

ABSTRACT  OF  #6 

In  this  paper,  the  performance  of  a  parallel  acquisition  scheme  for  a  direct  sequence 
spread  spectrum  receiver  operating  in  the  presence  of  data  is  presented.  The  effects  of  the 
data  modulation  on  key  quantities  such  as  probability  of  false  alarm,  probability  of  correct 
detection  and  mean-time  to  acquisition  are  analyzed  in  both  the  search  mode  and  the  lock 
mode,  and  a  search-lock  strategy  aimed  at  minimizing  the  acquisition  time  is  described. 

ABSTRACT  OF  #8 

The  acquisition  and  tracking  systems  of  a  spread-spectrum  receiver  are  probably  the 
most  critical  components  of  the  receiver,  since  if  they  fail  to  function  properly,  it  is  doubt¬ 
ful  that  the  desired  signal  can  be  successfully  detected.  This  means  that  the  affect  of 
interference  (such  as  jamming)  on  the  receiver  while  it  is  attempting  to  learn  the  correct 
phase  position  of  the  incoming  code  might  be  especially  harmful,  since  the  interference 
might  not  allow  the  receiver  to  acquire  the  signal.  To  address  this  problem,  a  narrow- 
band  interference  suppression  filter  is  used  to  enhance  the  performance  of  a  serial  search 
acquisition  scheme  for  a  direct-sequence  spread-spectrum  receiver.  Analytical  expressions 
for  the  probabilities  of  error  in  both  the  search  and  lock  modes  are  derived,  and  numerical 
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results  axe  used  to  illustrate  the  sensitivity  of  the  receiver  to  various  system  parameters.  It 
is  shown  that  the  presence  of  the  rejection  filter  can  significantly  improve  the  performance 
of  the  acquisition  system. 

5.  USE  OF  SPREAD  SPECTRUM  TECHNIQUES  IN  OPTICAL 
TRANSFORM  DOMAIN  PROCESSING 

The  results  are  discussed  in  paper  #7  in  Appendix  A.  This  paper  was  published  in 
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Figure  2.  Functional  Block  Diagram  of  the  System 
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Figure  5.  ?  versus  SNR,  showing  the  theoretical  and  exper¬ 
imental  results  for  notched  and  unnotched  single 
tone  interference  for  two  different  interference 


tones 


Figure  6 


dto/ir 


P  versus  the  notch  band-width  at  tvo  interference 
e 

tones  for  both  theoretical  and  experimental 
results . 
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Figure  8  Pfi  versus  SNR,  for  sweeping 
different  notch  band-widths. 


interference  and 
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ABSTRACT 

One  of  the  most  important  properties  of  direct-sequence  (DS)  spread 
spectrum  signals  is  their  low  probability  of  intercept.  Detection  of  the 
presence  of  such  signals,  as  well  as  characteristics  of  those  signals,  such 
as  center  frequency  and  chip-rate,  is  of  interest  to  an  intercept  receiver. 
However,  conventional  intercept  receivers  often  fail  when  they  operate  in  the 
presence  of  narrow-band  interference,  due  to  the  high  probability  of  false 
alarm. 

An  intercept  receiver  which  employs  transform  domain  processing  (TDP) 
is  constructed  at  a  center  frequency  of  200  Mhz  with  a  30  Mhz  bandwidth. 

This  receiver  detects  DS-BPSK  spread  spectrum  signals  in  the  presence  of 
narrow-band  interference  by  employing  adaptive  narrow-band  interference 
rejection  techniques.  The  improvement  in  the  system  performance  over  that  of 

* 
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conventional  detection  techniques  is  shown  by  presenting  experimental  meas¬ 
urements  of  probability  of  detection  versus  false  alarm  for  an  enhanced 
total  power  detector.  Also  presented  are  certain  results  corresponding  to 
double  frequency  line  detection  wherein  the  goal  is  not  just  signal  detection, 
but  also  carrier  frequency  estimation. 

The  receiver  employs  one  of  two  transform  domain  processing  techniques 
for  adaptive  narrow-band  interference  rejection.  In  the  first  technique  the 
narrow-band  interference  is  detected  and  excised  in  the  transform  domain  by 
using  an  adaptive  notch  filter.  In  the  second  technique,  the  interference  is 
suppressed  using  soft-limiting  in  the  transform  domain.  In  both  techniques 
the  real-time  Fourier  transformation  is  performed  using  surface  acoustic  wave 
chirp  filters  with  200  Mhz  center  frequency,  60  Mhz  bandwidth  and  60  us  delay. 
The  intercept  receiver  performance  is  studied  for  both  weighted  and  unweighted 
receiver  signals  at  various  inter ference-to-signal  and  signal-to-noise 
ratios,  and  in  the  presence  of  both  CW  and  finite-bandwidth  interference. 
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1.  INTRODUCTION 


Due  to  their  many  advantages,  spread  spectrum  communication  techniques 
have  become  attractive  to  both  the  military  and  civilian  communication 
sectors  [1-2].  Some  of  the  well-known  advantages  of  spread  spectrum  wave¬ 
forms  are  limited  immunity  to  narrow-band  interference,  transparency  of  the 
transmitting  signals  to  unintended  receivers  such  as  intercept  receivers, 
tolerance  to  multipath  effects,  and  code  division  multiple  access.  The 
capabilities  of  providing  such  important  properties  have  made  the  applica¬ 
tion  of  spread  spectrum  communication  especially  attractive  to  receivers 
operating  in  a  hostile  environment  where  low  probability  of  intercept  (LPI) , 
as  well  as  immunity  to  intentional  and  unintentional  interference,  is  highly 
desirable . 

Although  there  are  several  common  spread  spectrum  techniques,  such  as 
direct-sequence  (DS)  modulation,  frequency-hopping,  time-hopping  and  various 
combinations  of  the  above  waveforms,  the  only  technique  considered  in  this 
paper  is  DS,  generated  by  modulating  the  binary  information  with  a  wide-band 
pseudo-noise  (PN)  sequence,  and  further  up-converted  to  a  carrier  freqeuncy. 
In  what  follows,  the  application  of  real-time  transform  domain  processing 
(TDP)  [3]  for  adaptive  interference  suppression  in  an  intercept  receiver  is 
presented.  The  implemented  TDP  receiver  is  capable  of  adaptive  interference 
rejection  by  either  "excision"  or  "soft-limiting",  which  correspond,  respec¬ 
tively,  to  frequency  domain  notch  filtering  and  to  attenuating  appropriate 
interference  spectral  components  by  soft-limiting  [4-7].  In  the  results 
presented  below,  the  TDP  receiver  is  implemented  with  surface  acoustic  wave 
chirp  filters  having  60  Mhz  bandwidth,  60  us  delay,  and  a  center  frequency  of 


200  Mhz. 
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The  most  classical  technique  to  detect  the  presence  of  a  DS-BPSK 

signal  is  to  use  a  total  power  radiometer,  which  consists  of  a  BPF:  a  squaring 

device,  and  an  integrator.  The  detected  power  of  the  received  signal  is 

compared  to  an  external  threshold,  which  is  a  function  of  the  additive  thermal 

noise.  The  decision  on  the  presence  of  a  DS-SS  signal  is  made  if  the  measured 

detected  power  exceeds  the  threshold.  The  performance  of  such  a  detection 

technique  can  be  illustrated  by  presenting  the  probability  of  detection  (P.) 

d 

versus  the  probability  of  false  alarm  (Pr  ) ,  which  results  in  the  well-known 

ra 

receiver-operating-characteristic  of  the  detector. 

An  alternative  to  a  baseband  radiometer  is  a  feature  detection  scheme. 
The  feature  detection  techniques  are  signal  processing  configurations  capable 
of  detecting  particular  parameters,  or  features,  of  the  intercepted  spread 
spectrum  signal.  Feature  detectors  are  an  extremely  important  tool  in 
identifying  spread  spectrum  signals.  The  configuration  of  the  feature 
detector  depends  on  the  modulation  type  of  the  intercepted  signal  and  the 
desired  parameter  or  feature.  For  example,  if  one  replaces  the  integrator 
in  the  total  power  radiometer  with  a  narrowband  BPF  centered  at  twice  the 
frequency  and  follows  that  with  either  an  envelope  detector  or  a  square-law 
detector,  one  can  detect  the  so-called  "double-frequency  line." 

However,  the  operation  of  the  above  detection  technique  in  the  presence 
of  narrow-band  interference  typically  collapses  due  to  the  increase  in  the 
probability  of  false  alarm.  This  is  because  the  receiver  does  not  have  any 
apfu_o?u.  knowledge  regarding  either  the  presence  or  power  of  either  the  inter¬ 
cepted  signal  or  the  narrow-band  interference.  In  Section  2  of  this  paper,  we 
present  the  architecture  of  a  real-time  TDP  receiver  designed  to  alleviate 
the  above  problem,  and  Section  3  presents  experimental  results  of  the  enhanced 
total  power  radiometer  intercept  receiver.  In  Section  4,  the  use  of  the 
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adaptive  interference  suppression  scheme  is  shown  to  aid  in  the  determination 
of  the  center  frequency  of  the  intercepted  signal. 

The  results  of  Section  3  and  4  correspond  to  a  notch  filter  for  inter¬ 
ference  rejection.  Section  5  presents  the  application  of  transform  domain 

soft-limiting  for  interference  suppression.  Experimental  results  (P,  versus 

‘  d 

P^a)  in  the  presence  of  finite  bandwidth  interference  are  shown,  illustrating 
the  advantages  of  soft-limiting  over  the  notch-filter  technique.  Finally, 
Section  6  summarizes  the  results  of  using  the  TDP  intercept  receiver. 

2.  ARCHITECTURE  OF  THE  TRANSFORM  DOMAIN  PROCESSING  RECEIVER 

Real-time  realization  of  Fourier  transformation,  employing,  for 
example,  surface  acoustic  wave  (SAW)  chirp  filters  or  fast  Fourier  transform 
(FFT)  techniques,  has  opened  a  new  avenue  for  adaptive  signal  processing. 
However,  SAW  chirp  transformers  are  presently  a  far  superior  candidate  to  the 
FFT  counterpart  for  spread  spectrum  applications,  due  to  their  capability  for 
wideband  signal  processing  [8,9].  Hence,  in  what  follows,  SAW  filters  are 
exclusively  used. 

It  is  assumed  that  the  received  signal  is  a  composite  of  a  direct- 
sequence  binary-phase-shift-keying  (DS-BPSK)  signal,  narrow-band  interference, 
and  additive  thermal  noise.  Since  the  bandwidth  of  the  interference  is 
relatively  narrow  compared  to  the  intercepted  DS  spread  spectrum  signal, 
inspecting  the  transform  of  the  composite  signal  enables  the  processor  to 
decide  on  the  presence  and  center  frequency  of  the  interference.  Utilizing 
this  information,  two  techniques  are  implemented  for  rejecting  the  detected 
interference:  (1)  the  narrow-band  interference  is  adaptively  excised  by 

multiplying  the  transform  of  the  received  signal  with  an  appropriate  function 
containing  nulls  at  the  desired  locations;  (2)  the  interference  is  attenuated 
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to  the  desired  level  (such  as  noise  level)  employing  the  saturation  property 
of  an  amplifier. 

To  further  enhance  the  TDP  interference  rejection  schme,  the  received 
signal  is  time-weighted  employing  a  weighting  function  synthesizer.  This 
results  in  a  significant  reduction  of  the  interference  sidelobes.  The 
combination  of  interference  excision  and  time-weighting  attenuates  the  inter¬ 
ference  power  by  more  than  40  dB. 

A  general  block  diagram  of  the  implemented  intercept  receiver  is  given 
in  Figure  1.  The  receiver  performs  a  real-time  transformation  of  the 
received  signal.  This  is  followed  by  one  of  the  two  techniques  for  interfer¬ 
ence  rejection,  and  the  resulting  waveform  is  then  inverse  transformed. 
Finally,  the  output  of  the  inverse  transformer  is  used  as  the  input  to  the 
detection  device.  It  is  to  be  mentioned  that  to  measure  the  power,  a  digital 
integrator  using  a  6  bit  A/D  converter  with  a  12.5  Mhz  sampling  rate  was  used. 

Figures  2a  and  2b  present  the  schemes  for  real-time  transformation 

using  SAW  chirp  filters.  These  techniques  are,  respectively,  referred  to  as 

the  convolve-raultiply- convolve  (CMC)  and  the  multiply-convolve-multiply  (MCM) 

configurations.  The  impulse  response  of  a  SAW  chirp  filter  or  a  SAW  RAC 

2 

filter  (reflective  array  compressor)  is  given  as  cos(ujQt  i  Bt  ) ,  where 
and  B  are,  respectively,  the  center  frequency  and  the  slope  of  the  chirp 
signals.  Depending  on  the  sign  of  the  slope  (±  8)  the  chirp  signal  is  either 
an  up-chirp  or  a  down-chirp.  The  RAC  filters  employed  for  the  hardware 
implementation  of  the  receiver  are  all  down-chirp  RAC  filters. 

Figure  3  presents  the  transform  domain  processing  hardware  configura¬ 
tion.  For  simplicity,  consider  initially  just  the  desired  signal  in  the 
absence  of  noise  and  interference.  To  avoid  feed-through  of  the  input 


signal  given  by 
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S(t)  cos(ojQt),  0  1  c  1  T> 

to  the  RF  sections  of  the  TDP  receiver,  this  input  is  multiplied  by  a  600  Mhz 

carrier  frequency  prior  to  the  transform  domain  processing.  The  output  of 

the  modulator  (RF  mixer)  can  then  be  represented  by  S(t)  cos(4u>gt).  The 

desired  multiplying  chirp  component  employed  by  the  MCM  conf iguration  shown 

2 

in  Figure  2b  is  generated  by  allowing  the  chirp  signal  cos(a>Qt-6t  )  to 

modulate  a  400  Mhz  carrier.  Hence  the  received  signal,  now  centf -  I  at 

800  Mhz  is  modulated  by  the  desired  chirp  at  600  Mhz.  The  product  is 

2 

S(t)cos(a)Qt  +  8t  )  ,  and  this  is  used  as  the  input  to  a  RAC  filter  with  an 
impulse  response  given  by  cos((jj^t  -  Bt2). 

The  output  of  the  RAC  filter  is  given  by  [2] 


a0(t:) 


cos (u^t 


Bt2) 


S(t)  cos[2Btx]dT 


+  sinCu^t 


Bt2) 


S (t)  sin (2B  +  t) dt , 


(1) 


where  T  and  T^  are,  respectively,  the  duration  of  the  signal  and  the  delay 
of  the  RAC  filter  (duration  of  the  impulse  response).  Note  that  (1)  is  only 
valid  when  the  received  signal  is  fully  contained  within  the  RAC  filter 
(see  [2]).  The  output  of  the  RAC  is  fed  into  the  rejection  filter,  and 
this  is  followed  by  inverse  transformation  and  lo'J-pass  fitlering.  The 
output  of  the  inverse  transformation  is  translated  to  baseband  to  reduce  the 
hardware  complexity  of  the  detection  techniques. 

A  block  diagram  of  the  scheme  for  interference  excision  is  shown  in 
Figure  4.  The  transform  of  the  received  signal  is  split  into  two  paths; 
the  upper  path  is  modulated  for  interference  excision  by  employing  the  control 
signal  generated  by  the  lower  path;  the  lower  path  consists  of  an  envelope 
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detector  followed  by  logic  to  locate  the  position  of  the  interference,  as 
well  as  its  width.  Figure  5  presents  the  timing  of  the  upper  path's  trans¬ 
form  and  the  output  of  the  envelope  detector  for  three  consecutive  30  us 
periods.  (Note  that  30  us  is  taken  as  the  duration  of  each  input  segment 
to  be  processed.  This  is  because  the  RACs  have  time  delays  of  60  usee,  and, 
as  noted  above,  the  waveforms  to  be  transformed  must  be  fully  contained  in 
the  RAC  in  order  for  the  transform  to  be  valid.  Also,  as  discussed  in  [3j, 
the  time  delay  of  the  device  is  often  taken  to  be  twice  the  duration  of  the 
waveform  to  be  transformed.)  It  can  be  seen  that  the  envelope  detector 
introduces  a  constant  delay  on  the  signal’s  envelope;  therefore,  the  scheme 
has  to  compensate  for  this  delay  in  order  to  process  the  upper  path's  trans¬ 
form. 

The  first  option  for  compensating  the  delay  is  to  introduce  a  constant 
delay  in  the  upper  path.  This  would  require  applying  a  delay  line  with  a 
minimum  bandwidth  of  30  Mhz  operating  at  a  200  Mhz,  with  a  delay  of  about 
2  us.  A  good  candidate  is  a  surface  acoustic  wave  delay  line,  which  would 
meet  these  requirements;  however,  due  to  the  fixed  nature  of  such  a  delay 
line,  an  alternate  technique  is  implemented,  giving  the  receiver  the 
flexibility  to  have  a  variable  delay. 

In  the  second  technique,  it  is  assumed  that  the  interference  frequency 
does  not  change  within  the  two  consecutive  transform  periods,  or  faster 
than  7  Khz  (1/144  us).  Therefore,  the  decision  on  the  presence  and  location 
(frequency)  of  the  detected  interference  in  period  it  1  can  be  applied  to  the 
interference  within  period  it 2.  Similarly,  the  information  obtained  from 
period  it 2  can  be  used  to  process  period  it 3,  etc. 

Figure  6  presents  a  detailed  block  diagram  of  the  implemented  excision 
scheme.  The  output  of  the  envelope  detector  is  followed  by  a  hard  decision 
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scheme  which  incorporates  a  threshold.  The  30ps  output  of  the  envelope 
detector  is  sampled  at  a  rate  of  25  Mhz,  where  the  sampling  region  is 
controlled  by  the  read/write  (R/W)  signal  initiated  from  the  system  timing. 

Each  sampled  value  is  compared  to  the  pre-assigned  threshold  and  a  hard 
decision  is  made  (0  or  1,  depending  on  whether  the  threshold  has  been  exceeded) . 
Each  sampled  point  is  assigned  an  address  starting  from  (0000,0000,000)  to 
(1111,1111,111);  therefore,  each  address  value  will  correspond  to 
a  frequency  component.  The  information  regarding  the  first  transform  period 
(period  Tl;  Figure  5a)  is  stored  at  each  point's  assigned  address  in 

memory  it  A.  At  the  end  of  period  Tl,  memory  if  A  will  have  the  information 

regarding  the  presence  of  interference  in  period  Tl.  The  next  period 

(period  T2)  is  sampled  at  25  Mhz  with  the  same  addressing  assignments; 

however,  this  information  is  stored  in  memory  if B. 

While  storing  the  information  regarding  period  T2,  information  stored 
in  memory  it  A  (corresponding  to  period  Tl)  is  retrieved,  and  the  initial  read¬ 
ing  address  is  advanced  from  (0000,0000,000)  to  (abed, efgh, 000)  where  abed  and 
efgh  are  adjustable.  Therefore  the  timing  of  the  retrieved  information  is 
advanced  with  respect  to  period  T2  (Figure  5a)  and  the  delay,  due  to  the  envelope 
detector  and  the  electronics,  is  compensated.  While  observing  period  T3,  decis¬ 
ions  are  stored  in  memory  if  A  and  information  stored  in  memory  It  B  is  retrieved 
using  the  same  initial  address.  The  retrieved  information  from  the  memories  is 
shown  in  Figure  5b,  which  also  indicates  the  location  of  interference  within 
the  30  us  transform  period.  This  information  (narrow  pulses)  is  used  to 
generate  gating  pulses  with  an  adjustable  width  (excision  pulses) .  These 
pulses  overlap  the  interference  in  the  transform  period,  as  can  be  seen 
from  Figure  5b.  It  is  necessary  to  adjust  the  initial  address  for  retrieving 
the  information,  so  that  the  excision  width  overlaps  the  region  of  the 


8 


interference's  transform  to  be  excised.  Employing  this  scheme,  the  receiver 
das  complete  flexibility  to  align  the  excision  pulses  with  the  upper  path's 
interference.  However,  this  would  not  be  possible  if  a  constant  delay  was 
introduced  in  the  upper  path. 

The  excision  pulse  is  used  to  control  the  RF  switch  in  the  upper  path, 
therefore  excising  the  interference.  The  amount  of  excision  depends  on  the 
interference  bandwidth.  In  the  implemented  scheme,  the  excision  width  has 
to  be  manually  adjusted.  Figure  7a  presents  the  output  of  the  transformer 
for  a  composite  of  5  Mhz  chip-rate  DS-BPSK  signal  and  a  pulsed  CW  interfer¬ 
ence  with  a  1.625  Mhz  pulse-rate.  Figure  7b  shows  the  same  transform  after 
interference  excision.  Figures  8a  and  8b  present  the  spectrum  of  the  inverse 
transformer's  output  before  and  after  interference  excision,  respectively.  It 
can  be  seen  that  the  components  exceeding  the  threshold  have  been  excised. 

Figures  9a  and  9b  present  the  output  of  the  transformer  for  a  single 
CW  interference  before  and  after  interference  excision.  Figures  10a  and  10b 
present  the  corresponding  time  domain  signals  at  the  output  of  the  inverse 
transformer;  the  interference's  residue  due  to  the  side-lobes  can  be  seen  in 
Figure  10b.  To  optimize  the  receiver's  performance,  it  is  imperative  to 
minimize  the  interference's  power  as  much  as  possible,  with  a  minimum  loss 
of  the  desired  signal.  To  reduce  the  interference's  residue,  two  options 
are  possible:  (1)  increase  the  excision  width,  which  would  not  always  be 
feasible  in  the  event  of  multiple  interference  or,  (2)  time-weight  the 
received  signal. 

To  reduce  the  contribution  of  the  interference  sidelobes  employing 
the  weighting  technique,  either  the  received  signal  or  the  RAC  filter  must  be 
weighted.  Since  weighting  of  a  RAC  filter  is  both  an  expensive  and  a  permanent 
process,  the  input  signal  is  time-weighted  instead.  The  weighting  function 
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that  is  used  is  the  6  terra  Blackman-Harris  window  [10]  given  by 

W(n)  =■  aQ  -  cos  n]  +  a2  cos  [  2n]  -a3  cos  [^-  3n]  ,  (2) 

n  -  0,1, _ (N-l) , 

where 

aQ  =  0.35875, 
a^  =  0.48829, 
a2  =  0.14128, 

a3  -  0.01168,  (3) 

and  N  is  the  number  of  samples  within  the  signal  period. 

The  30  ps  long  weighting  function  is  sampled  at  25  Mhz,  resulting  in 
750  samples.  The  analytical  sample  values  are  digitized  and  the  8  bit 
numbers  are  stored  in  an  EPROM.  Figure  11  presents  the  implemented  scheme 
for  the  reproduction  of  the  weighting  function.  To  synthesize  the  stored 
values,  the  contents  of  the  EPROM  are  transferred  to  a  high  speed  RAM 
(since  a  commercial  EPROM  cannot  be  operated  at  25  Mhz) .  The  8  bit  sample 
values  are  retrieved  from  the  RAM  at  a  rate  of  25  Mhz,  and  converted  to  the 
ECL  level.  This  is  followed  by  a  high  speed  ECL  analog-to-digital  (A/D) 
convertor.  Finally,  the  output  of  the  A/D  convertor  amplitude  modulates 
the  received  signal  employing  an  RF  mixer.  The  timing  signals  originating 
at  the  TDP  receiver  control  the  beginning  of  the  30  ps  weighting  function 
period,  so  that  it  overlaps  the  TDP's  multiplying  chirp. 

The  interference  excision  technique  uses  the  TDP  receiver  to  excise  a 
desired  interference  bandwidth.  However,  in  notching  the  interference 
components,  a  part  of  the  signal  overlapping  the  same  frequency  and  band¬ 
width  is  also  excised.  The  optimum  excision  width  is  a  function  of  the 
interference  frequency  and  bandwidth.  Therefore  the  interference  band¬ 
width  increases,  so  does  the  excision  width  and  the  excised  signal  power. 

An  alternative  to  this  is  to  employ  a  soft-limiter.  Through  soft-limiting, 
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the  interference  components  can  be  attenuated  to  a  desired  level  in  the 
transform  domain,  thus  resulting  in  the  signal  components  having  similar 
distortion. 

Two  possible  soft-limiting  techniques  are  presented  here.  The  first 
technique  is  shown  in  Figure  12a,  where  the  transform  of  the  received  signal 
is  split  into  two  paths;  the  lower  path  is  followed  by  an  envelope  detector 
and  further  sampled  using  a  high  speed  6  bit  analog-to-digital  (A/D)  con¬ 
vertor.  The  6  bit  values  correspond  to  the  spectral  magnitudes  of  the 
transform.  This  information  is  processed  and  a  signal  is  generated  and  used 
to  amplitude  modulate  the  transform  in  the  upper  path.  This  technique  has 
to  compensate  for  the  delay  introduced  by  the  envelope  detector  and  the 
electronics;  this  compensation  takes  the  same  form  as  the  one  discussed 
above  for  the  excision  technique. 

Figure  12b  presents  the  second  technique  for  soft-limiting.  In  this 
scheme,  the  transform  of  the  received  signal  is  followed  by  a  variable 
attenuator  and  a  high-gain  amplifier.  Since  the  interference  is  narrow-band, 
but  high  level,  it  is  clipped  due  to  the  saturation  of  the  amplifier,  and 
the  level  of  clipping  is  controlled  by  the  front-end  variable  attenuator. 

This  amplifier  is  followed  by  a  second  attenuator,  used  to  adjust  the  power 
level  to  the  inverse  transformation  stage. 

3.  TOTAL  POWER  RADIOMETER 

A  general  block  diagram  of  the  measurement  setup  is  shown  in  Fig¬ 
ure  13.  The  received  signal,  a  composite  of  a  direct-sequence  BPSK  signal, 
additive  thermal  noise,  and  narrow-band  interference,  is  first  gated  in 
time.  The  time-gating  employs  either  a  30  us  long  rectangular  weighting 
function,  or  the  6  term  Blackman  Harris  window  function,  analytically 
given  by  Equation  (2). 
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As  can  be  seen  from  Figure  13,  the  5  Mhz  chip  rate  DS-BPSK  signal 
is  generated  through  a  spread-spectrum  signal  generator,  and  the  interfer¬ 
ence  is  either  a  sinewave  or  a  DS-BPSK  signal  with  a  variable  (low)  chip-rate, 
the  bandwidth  of  which  is  taken  to  be  its  main  lobe.  (Figure  3  does  not  show 
the  pure  sinewave  interferer.)  The  results  to  be  presented  in  both  this 
section  and  Section  4  correspond  to  the  CW  sinewave  interference,  whereas 
the  results  in  Section  5  correspond  to  a  finite  bandwidth  interferer. 

To  generate  the  additive  thermal  noise,  a  wide-band  noise  source  is  followed 
by  a  bandpass  filter  with  a  60  Mhz  bandwidth,  where  60  Mhz  is  the  front-end 
bandwidth  of  the  transform  domain  processing  receiver.  The  respective 
powers  of  signal,  noise  and  interference  are  measured  following  the  time¬ 
gating  of  the  composite  signal.  Using  two  step  attenuators,  the  signal-to- 
noise  (SNR)  and  the  interf erence-to-signal  ratios  (I/S)  can  be  appropriately 
adjusted. 

The  intercepted  DS-BPSK  is  assumed  to  have  a  3  Mhz  chip  rate,  resulting 
in  a  10  Mhz  mainlobe  centered  at  -  209  Mhz.  The  interference  components 
are  also  assumed  to  be  within  the  bandwidth  of  the  DS-BPSK  signal.  The 
composite  received  signal  is  continuous;  however,  due  to  the  nature  of  TDP 
(processing  a  30  us  period  every  144  us),  the  output  of  the  receiver  will 
also  be  a  30  us  signal  every  144  us.  The  output  of  the  TDP  receiver  is 
followed  by  the  total  power  radiometer.  The  integrator  samples  the  output 
of  the  detector  at  a  rate  of  12.5  Mhz,  and  accumulates  the  6  bit  values 
over  the  30  us  period.  At  the  end  of  the  integration  time,  the  detected 
power  (18  bit  value)  is  compared  to  an  external  threshold.  The  number  of 
positive  decisions  is  recorded  within  10^  (30  us  signal  period)  observations. 
The  performance  of  the  radiometer  is  shown  by  presenting  the  probability  of 
detection  and  false  alarm  given  as 
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P,  =  n,/M  and  Pp  =  nc  /M  , 
d  d  fa  fa 

where  M  is  the  number  of  observations  (10^)  and  n^  amd  m^  are,  respectively, 
the  number  of  positive  detections  in  the  presence  and  in  the  absence  of  the 


DS-BPSK  signal. 

The  performance  of  the  receiver  utilizing  the  adaptive  interference 

excision  scheme  is  shown  in  Figure  14,  which  presents  the  measured  P,  and  Pc 

d  f  a 

for  an  interference  to-signal  ratio  (I/S)  of  24  dB  with  the  excision  filter 
on,  and  I/S  ratios  of  4  and  5  dB  with  the  excision  filter  off.  The  inter¬ 
ference  in  this  case  is  a  fixed  frequency  sinewave.  Figure  15  shows  results 
for  a  swept  frequency  sinewave,  where  the  interferer  is  swept  over  the 
10  Mhz  signal  bandwidth  at  a  I/S  =  24  dB  (excision  scheme  on)  and  I/S  =  2  and 
3  dB  (excision  scheme  off) .  It  can  be  seen  from  Figures  13  and  14  that  in 
both  cases,  the  interference  excision  scheme  has  enhanced  the  receiver's 


performance  by  an  excess  of  20  dB. 

Since  the  TDP  excises  the  main  lobe  and  a  few  additional  sidelobes  of 


the  interference,  it  is  anticipated  that  if  the  received  signal  is  time- 
weighted  prior  to  the  transform  domain  processing,  further  improvement  in 
the  receiver's  performance  will  result.  In  particular,  it  is  expected  that 
at  high  I/S  ratios  time-weighting  would  be  extremely  advantageous.  Figure 
16  presents  the  performance  of  the  receiver  for  weighted  and  unweighted 
received  signals  at  SNR  =  0  dB  and  I/S  ratios  of  28  and  33  dB.  It  is  seen 
that  the  system  using  the  weighted  signal  outperforms  the  one  employing 
rectangular  weighting. 


DOUBLE  FREQUENCY  LINE  DETECTION 


Figure  17  presents  a  block  diagram  of  the  experimental  setup  for 
the  double  frequency  line  detection.  The  composite  signal  is  generated 


13 


by  combining  the  DS-BPSK  signal,  additive  thermal-noise  and  narrow-band 
interference  as  before.  Since,  ordinarily,  one  does  not  know  the  center 
frequency  of  the  received  signal,  a  BPF  with  a  variable  center  frequency  is 
needed  following  the  square-law  device.  However,  to  simplify  the  experimental 
setup,  a  fixed  bandpass  filter  with  a  500  Khz  bandwidth  is  employed. 

Figures  18-20  illustrate  the  performance  of  the  double  frequency  line 
detection.  These  curves  are  somewhat  different  from  the  curves  of  Section  4, 
in  that  Pp  is  no  longer  plotted  against  P^  ,  hut  rather  is  plotted  against 
the  probability  of  detecting  the  signal  resulting  from  the  crossproduct  of 
the  transmitted  waveform  and  the  interference.  The  reasons  for  this  are  as 
follows:  For  the  parameter  values  used  in  the  measurements,  the  P  was 

I  3 

very  small  when  the  TDP  rejection  scheme  was  used.  Therefore,  it  was  decided 
to  emphasize  the  probability  of  correctly  identifying  the  actual  double 
frequency  of  the  intercepted  signal  (to  within  the  bandwidth  of  the  narrow- 
band  BPF)  relative  to  the  probability  of  falsely  identif yi»-0  the  sume 
frequency  generated  by  the  crossproduct  of  the  signal  and  interference  at 
the  output  of  the  squaring  device.  Note  that,  in  principle,  there  is  also 
present  at  the  output  of  the  squarer  a  term  located  at  the  double  frequency 
of  the  interference.  However,  when  the  TDP  system  was  employed  in  the 
receiver,  this  latter  term  was  found  to  be  significantly  smaller  than  the 
cross-term  referred  to  above. 

Figure  18  shows  the  performance  of  this  system  for  two  different 
sinewave  interferers,  one  located  at  211.65  Mhz,  and  the  other  located  at 
212.65  Mhz.  Since  the  center-frequency  of  the  DS-BPSK  signal  is  209.65  Mhz, 
it  is  seen  that  the  system  experiences  greater  degradation  when  the  inter¬ 
ference  is  closer  to  the  desired  signal. 


14 

Figure  19  presents  the  performance  of  the  detector,  both  with  and 
without  interference  excision,  for  I/S  =  1  dB  (excision  off),  and  I/S  =  20  dB 
(excision  on).  It  can  be  seen  that  the  TDP  excision  scheme  has  again  improved 
the  receiver's  performance  by  an  excess  of  20  dB.  Also  shown  on  this  figure 
is  the  relative  performance  of  the  system  when  detecting  the  double  frequency 
line  of  the  interference  when  TDP  is  not  employed  with  I/S  «  20  dB. 

Figure  20  presents  similar  results  when  the  detector  operates  in  the 
presence  of  pulsed  narrow-band  interference  with  a  1.617  Mhz  pulse  rate. 

For  this  latter  case,  there  are  three  interference  components  within  the  main 
lobe  of  the  DS-BPSK  signal.  The  abscissa  now  corresponds  to  the  probability 
of  detecting  the  crossproduct  between  the  DS-BPSK  signal  and  the  main  spectral 
component  of  the  interference.  The  measured  experimental  performance  corres¬ 
ponds  to  I/S  =  10  dB,  and  the  center  frequencies  of  signal  and  interference 
are  the  same  as  for  the  results  presented  in  Figure  18.  It  can  be  seen  that 
without  the  adaptive  interference  excision,  the  performance  is  significantly 
degraded,  and  for  any  given  P^,  the  probability  of  incorrectly  identifying 
the  cross  term  is  approximately  unity. 

4.  FINITE-BANDWIDTH  INTERFERENCE  REJECTION 


Section  3  illustrated  the  performance  of  the  intercept  receiver 
operating  in  the  presence  of  CW  jammers,  and  employing  TDP  interference 
excision.  In  this  section,  the  performance  of  the  intercept  receiver  operating 
in  the  presence  of  finite  bandwidth  interference  is  presented  for  each  of  two 
rejection  techniques:  (1)  interference  excision  and  (2)  soft-limiting,  both 
employing  transform  processing. 

The  experimental  setup  is  that  shown  in  Figure  1.  The  composite 
signal  is  generated  by  combining  a  5  Mhz  chip-rate  DS-BPSK  signal  with  the 
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finite  bandwidth  interference  and  thermal  noise.  The  PN  sequence  used  for 

32 

generating  the  interference  has  a  period  of  2  -1.  The  soft-limiting  tech¬ 

nique  attenuates  the  interference's  spectral  components  to  a  desired  level. 

In  the  scheme  implemented  here,  this  is  achieved  by  employing  the  saturation 
property  of  an  amplifier.  Saturation  of  the  amplifier  is  a  non-linear 
operation;  however,  at  high  inter ference-to-signal  ratios,  the  signal 
component  is  not  significantly  affected,  sicne  the  saturation  is  in  transform 
domain  and  the  signal's  spectral  component  is  much  lower  than  the  inter¬ 
ference's  component.  The  performance  of  the  receiver  employing  the  soft 
limiting  technique  is  directly  proportional  to  the  amplifier's  saturation 
level,  the  interference's  power  and  bandwidth,  and  the  signal-to-noise 
ratio.  Prior  to  any  measurements,  the  operation  of  the  receiver  is  optimized 
by  adjusting  the  two  variable  attenuators  shown  in  Figure  13. 

Figures  21  and  22  present  the  measured  performance  of  the  excision 
technique  and  the  soft  limiting  scheme,  respectively,  for  various  interfer¬ 
ence  bandwidths.  Figure  23  shows  a  direct  comparison  of  the  two  schemes,  and 
the  advantages  of  soft-limiting  over  the  excision  technique  can  be  seen.  The 
main  complications  with  finite  bandwidth  interference  is  that  the  excision 
bandwidth  cannot  be  increased  without  losing  a  significant  amount  of  inter¬ 
cepted  signal's  power,  especially  when  operating  in  the  presence  of  multiple 
interference . 

5.  SUMMARY  AND  CONCLUSION'S 

This  paper  has  demonstrated  the  advantages  of  using  transform  domain 
processing  when  trying  to  intercept  a  DS-BPSK  signal  in  the  presence  of 
strong  narrow-band  interference  and  additive  thermal  noise.  The  TDP 
receiver  presented  here  is  implemented  at  a  center  frequency  of  200  Mhz  with 
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a  30  Mhz  processing  bandwidth.  The  receiver  is  experimentally  analyzed  and 
its  performance  measured  under  various  signal-to-noise  ratios  in  the 
presence  of  both  CW  tone  and  finite  bandwidth  interference.  The  performance 
of  the  receiver  is  further  enhanced  by  employing  time-weighting  at  the  front 
end  of  the  receiver;  this  is  especially  advantageous  at  high  I/S  ratios, 
where  the  side  lobes  have  a  strong  contribution. 

The  two  interference  rejection  techniques  that  have  been  implemented 
and  experimentally  analyzed  are  interference  excision  and  soft  limiting. 

Both  techniques  were  shown  to  significantly  improve  the  performance;  however 
the  measured  results  indicate  tnat  in  the  presence  of  finite  bandwidth 
interferers,  the  soft-limiting  technique  outperforms  the  excision  technique. 

The  performance  of  the  receiver  can  be  further  improved  by  optimizing 
the  excision  width,  the  weighting  function  and  the  level  of  the  soft-limiting 
based  on  the  interference  pwoer  and  center  frequency. 
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Figure  1  General  block  diagram  of  the  implemented 
receiver . 
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Figure  2a.  Convolve-multiply-convolve  real  time 
transformation  technique. 


Figure  2b.  Multiply-convolve-multiply  real  time 
transformation  technique. 


Real  Time  Transformation  Techniques  Employing 
Chirp  Filters 


Figure  3.  Transform  domain  processing  hardware  configuration. 


Transformer  Output 
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Figure  7a.  Composite  transform  of  a  5  Mhz  chip  rate  DS-BPSK 
and  a  pulsed  CW  interference  with  1.625  Mhz  pulse 
rate . 


Figure  7b.  Composite  transform  of  a  5  Mhz  chip  rate  DS-BPSK 
and  a  pulsed  CW  interference  with  1.625  Mhz  pulse 
rate  following  interference  excision. 
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Figure  10a.  30  us  CW  input  signal. 


Figure  10b.  Output  of  the  receiver  for  a  30  us  input  given  by 
2.22a  using  excision  technique. 
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Timing  for  the  Weighting 


Figure  11.  General  block  diagram  of  the  implemented 
weighting  function  synthesizer. 


Figure  12a.  Adaptive  soft-limiting  technique  employing  TDP 


un 


Figure  12b.  Implemented  TOP  sof t-- L imit ing  using  amplifiers. 
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Figure  14.  Measured  performance  for  I/S  =  4  dB  and  5  dB,  excision 
off  and  I/S  =  24  dB,  excision  on. 
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Measured  performance  at  I/S  =  31  dB  and  SNR 
for  excision  and  soft-limiting  technique  at 
interference  bandwidth  of  240  Khz. 
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Abstract 

Because  a  compressive  receiver  can  accomplish  signal  detection  in  a 
manner  analogous  to  that  of  a  parallel  bank  of  narrowband  filters,  it  has 
obvious  applications  in  the  area  of  signal  interception.  That  is,  it  can 
be  used  not  only  to  perform  the  most  basic  function  of  signal  detection, 
but  can  also  be  used  to  provide  information  on  frequency  location. 
However,  it  must  be  determined  if  there  is  any  significant  loss  in 
detection  performance  that  accompanies  this  added  feature  of  frequency 
estimation. 

This  paper  will  attempt  to  answer  that  question  by  comparing  the 
perf ormance  of  a  compressive  receiver  to  that  of  wideband  radiometer.  The 
waveforms  to  be  detected  are  either  sinewaves  or  phase-shift  keyed  signals, 
where  the  latter  waveforms  can  represent  either  a  narrowband  binary  PSK 
signals  or  direct  sequence  spread  spectrum  waveforms. 


This  work  was  partially  supported  by  the  US  Army  Research  Office  under 
Grant  DAA  L03-86-K-0092. 


I.  INTRODUCTION 

As  is  well  known,  there  are  many  approaches  available  for  the 
interception  of  digital  waveforms  [l]-[7],  There  are  also  a  variety  of 
functions  of  an  interception  system,  ranging  from  the  most  basic  one  of 
signal  detection,  to  more  sophisticated  ones  involving  feature  extraction. 
Detection  is  often  accomplished  by  either  a  wideband  radiometer  or  a 
channelized  radiometer,  where  the  latter  receiver  has  the  obvious  advantage 
over  the  former  receiver  of  providing  information  on  frequency  location  in 
addition  to  making  a  decision  on  whether-or-not  signal  is  present.  As  a 
practical  means  of  implementing  a  channelized  filter  bank,  a  compressive 
receiver  can  be  used  [7 1 . 

However,  if  a  compressive  receiver  is  to  be  considered  in  place  of 
a  radiometer,  it  must  be  determined  how  the  detectability  of  one  compares 
to  the  detectability  of  the  other.  That  is,  it  must  be  determined  if  the 
extra  frequency  resolution  of  the  compressive  receiver  results  in  decreased 
detectability. 

In  this  paper,  we  attempt  to  answei  the  question  posed  above. 
Performance  analysis,  which  consists  of  derivations  of  the  probability  of 
false  alarm  and  the  probability  of  detection  for  both  a  tone  signal  and  a 
binary  PSK  signal,  is  presented  for  both  systems,  and  comparisons  between 
the  two  are  made.  Also,  while  the  results  are  presented  for  a  binary  PSK 
waveform,  it  will  be  clear  that  the  same  analysis  applies  for  a  direct 
sequence  spread  spectrum  signal. 

In  ..he  next  section,  the  system  model  of  a  compressive  receiver  is 
reviewed  and  the  expressions  for  the  probability  of  false  alarm  and  the 
probal  lity  of  detection  are  derived.  Section  III  presents  the 
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corresponding  expressions  for  the  radiometer.  Numerical  results  for  the 
performance  of  both  receivers  are  provided  in  Section  IV,  and  conclusions 
are  presented  in  Section  V. 


II.  ANALYSIS  OF  COMPRESSIVE  RECEIVER 

A  block  diag -am  of  the  compressive  receiver  is  shown  in  Fig.l.  The 
input  X( t )  to  the  compressive  receiver  is  described  by  one  of  the  following 
two  hypotheses: 

H  :  X(t)  =  n  (t) 

0  w 

or 

H,  :  X(t)  =  fit)  +  n  (t)  , 

1  w 

where  fit)  is  the  signal  waveform  and  n  (t)  is  additive  white  Gaussian 
noise  of  two-sided  spectral  density  Vq/ 2.  This  input  is  multiplied  by  a 
periodic  scanning  waveform  ait),  and  the  product  sit)  is  used  as  an  input 
to  a  time-invariant  filter  hit).  Over  one  scan  period,  ait)  is  represented 


ait)  =  cos 


V 


-  Pt 


PT(r) 


(1) 


where  the  puise  function  P  (x )  is  defined  as 

a 


P  (x)  = 
a 


1 

0 


,  0  i  x  s  a 
,  elsewhere  . 


(2) 


The  impulse  response  hit)  of  the  filter  is  a  positive-going  chirp  for  T} 
seconds,  i.e., 


hit)  =  cos |  u  t  +  £t 


P  it)  . 
1 


(3) 
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In  the  following  analysis,  we  assume  that  T. 

We  want  to  obtain  the  receiver  output  in  one  scan  perioH  T.  The 
input  s(t)  is 

s(t)  =  X(t)  cos ^  «Qt  -  /3t2  j  PT(t)  ,  (4) 

and  the  output  g(t)  is 

r  ” 

g(t)  =  s(x)  h(t  -  t)  dx  .  (5) 

J  —00 

From  (2)-(5),  it  is  easily  shown  that,  for  any  T  ,  T  and  t, 

mlnft.T) 

g(t)  s  i-cos(wQt  +  0t2)  J  X(x)  cos(2gtx)  dx 

max  ( 0,t  -  T  ) 

1 

mlnft.T) 

+  i-sin(u>0t  +  0t2)  j"  X(x)  sin(20tx)  dx  ,  (6) 

mi  x(0,t-T  ) 

1 

where  the  approximation  results  from  ignoring  the  double  frequency  terms 
during  the  integration. 

Assume  that  the  output  git)  of  the  compressive  receiver  is  fed  into 
an  envelope  detector,  and  the  output  e(t)  of  the  envelope  detector  is 
sampled  (perhaps  multiple  times)  in  order  to  make  a  decision  as  to  whether 
or  not  the  signal  fit)  is  present.  This  is  illustrated  in  Fig.  2. 

A.  Tone  Signal 

We  are  initially  interested  in  the  case  where  fit)  is  a  pure  tone, 

i . e. , 

fit)  =  A  co$|  Wjt  +  0  j  ,  (7) 
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where  A  is  a  constant  amplitude  and  0  is  a  random  phase.  When  the  signal 
is  present,  (6)  gives 


git)  £  ^  cos(«oQt  +  3t2) 


sin 


Upt-o^H-e 


-  s  1 


23t 


n  (2/31-0^  Ht-T^-e 


u. 


^  siniu  t  +  3t2) 
4  0 


cos 


i^t-u^t-B 


-  cos 


i^t-tJ^it-TjJ-e 


23t  -  o) 


+  n(t) 


(8) 


The  noise  component  nit)  is  given  by 


nit)  =  lr  n  it)  cosiu  t  +  3t2)  +  4  n  (t)  siniu  t  +  3t2)  , 
2  c  0  2  s  0 


(9) 


where 


n  it)  = 
c 


t-T 


n  (t)  cos(23tx)  dr 
w 


(10) 


and 


n  it)  = 
s 


t-T 


n  ix)  sin(23tx)  dx 
w 


ill) 


Both  n^i t )  and  n^it)  are  zero-mean  Gaussian  random  processes.  If  it  is 

assumed  that  23tT,»27r,  then  the  variance  of  either  4  n  it)  and  4  n  it) 
1  2  c  2  s 

is  approximately 


%  =  V/16 • 


(12) 


and  n^it)  and  n^i t )  are  statistically  independent  of  one  another. 


The  output  of  the  envelope  detector  due  to  fit)  is  given  by 
AT, 


e  ( t )  = 


1 


sin [ (23t  -  w  )T  /2] 
( 23 1  -  ^lTj/2 


(13) 
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With  this  model,  an  output  sample  of  the  envelope  detector  will  have  a 
Rayleigh  density  if  the  signal  is  absent,  and  will  have  a  Rician  density  if 
the  signal  is  present.  Assuming  the  final  decision  is  made  on  the  basis  on 
a  single  sample,  the  probabilities  of  false  alarm  and  correct  decision  are 
given  by 


n 


(14) 


and 


Pd  =  Q 

respectively,  where  y  is  threshold  of  the  test,  t  is  the  sampling  time, 
and  Q(a,  b)  is  defined  as 


jr _ 

’  <r 


(15) 


Q(a,  b)  * 


(2  2  \ 

-  -  2~X~  j  I0(ax)  dX  ' 


(16) 


In  (16),  IQ(x)  is  the  zeroth-order  modified  Bessel  function  of  the  first 


kind. 

Rather  than  use  a  single  sample,  which  would  be  appropriate  if  the 
frequency  of  the  tone  was  known,  multiple  samples  of  the  output  are  needed 
to  make  the  final  decision  when  the  carrier  frequency  is  unknown.  To  be 
specific,  assume  that  N  independent  samples  of  the  envelope  detector  output 
can  be  accumulated  in  T  -  T  seconds.  A  decision  to  detect  the  signal  as 
being  present  is  made  if  any  one  of  the  N  samples  exceeds  the  threshold. 
Then  the  overall  probability  of  false  alarm  is  given  by 

P_  =  1  -  (1  -  P.  )N  ,  (17) 

FA  fa 


where  P. 

f  a 


is  given  by  (14),  and  the  overall  probability  of  correct  decision 
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is  given  by 


P 


D 


-  n  “ -  V  • 

l  =1  J 


(18) 


In  (18),  P  .  is  given  by  (15)  evaluated  at  the  sampling  time  t  .  Also,  in 

a.  s. 

l  1 

(17),  it  is  assumed  that  signal  is  absent  on  all  N  samples,  while  in  (18) 

it  is  assumed  that  signal  is  present  on  all  N  samples.  Note,  however,  that 

if  e,(t  )  =  0,  then  (15)  reduces  to  (14),  and  hence  (18)  also  applies  if 

1  s 

only  some,  but  not  all,  of  the  N  samples  contain  signal. 

To  determine  N,  consider  again  the  noise  at  the  output  of  the 

compressive  receiver.  This  can  be  expressed  as  in  (9),  or,  equivalently, 
as 


'.(t)  = 


1  n  (t-x)  cosLjt-ti  -  £(t-T)Z  cos  u  x  +  |3x2  dx 

w  [0  J  L  °  I 


1 

2 


1  f- 

n  (t-x)  cos  urtt  - 
w  0 

1 


0t  +  2£tx  dx 


(19) 


Hence,  the  autocorrelation  function  of  n(t)  is  given  by 


Vl 

R  (x)  =  -  ~  R.(x)  cos 
n  16  1 


jwjrl  +  0t2  +  ^j(t)J  , 


where 


f 

sinOxfT^-frf)] 

Rl(x)  -  | 

0xTi 

0 


M  *  T 

»  1 

,  elsewhere 


(20) 


(21) 


and 


7 


#j{t)  s  -tan 


9  ' 

cos(2ptTi)  -  cos(2pt  ) 
sin(2PMTi)  -  sin(2ptZ) 


(22) 


We  see  that  R  (t)  =0  if 
n 

PtCTj  -  t)  =  nir  ,  (23) 

where  n  is  any  integer.  Solving  (23)  for  T  yields  (approximately)  the  two 
solutions 


_  nrr  n  n 
T  1  PT  *  pT  * 

2 

where  the  approximation  holds  if  pT^  »  4nn.  From 
samples  of  g(t)  uniformly  spaced  at  intervals  of  rc/pT^ 
independent.  Hence,  in  T  -  seconds,  approximately 


(24) 

(24),  notice  that 
are  approximately 


N  =  pT  (T  -  T  )/n  =  2BTj  (25) 

independent  samples  can  be  accumulated,  where  2B  is  the  bandwidth  (in  Hz) 
to  be  searched  by  the  compressive  receiver. 


B.  Binary  PSK  Signal 

Suppose  that  f(t)  is  a  binary  PSK  signal  given  by 


M-  i 


f(t)  =  A 


I 

1  =o 


diPT, 


(t  - 


iT,  )  cos(w,t 
b  1 


+  0) 


(26) 


where  the  notation  is  the  same  as  that  of  the  previous  section  except  that 
is  the  duration  of  a  data  bit,  and  d.  is  equally  likely  to  be  ±1.  Also, 
we  let  MT^=T  and  mT  =T  ,  where  K1  and  m  are  integers.  If  we  again  assume 
that  Tj  ^  t  i  T  and  ti.e  signal  f(t)  is  present,  it  can  be  shown  that, 
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analogous  to  (8),  g(t)  is  given  by 


g(t)  =  ^  gc(t)  cos(uQt  +  /3t2)  +  -  gs(t)  sin(uQt  +  0t2)  +  n(t)  , 


where  n(t)  is  again  given  by  (9), 


.  (  sin[(2jSt-u  )(j-m)T  -  0)  -  sin[(23t-o>  )(t-mT,  )  -  0] 

A  I  ,  lb  lb 

~  2  \  j-m-1  2gt  -  u 


j-2  sint  (23t-u,  )(i+l  )T,  -  0]  -  sin[ (2£t-w,  )iT,  -  0] 

E,  lb  lb 

i  2/3t  -  u 

1  =  J-m  M  1 


+  dM 


sin [ ( 2/3t— ) t  -  0]  -  sin[(23t-Uj)( j-DT^  -  0] 
2p t  -  «, 


gs(t)  =  -  ? 


{  d j-m-1 


cost  (23t-u, )( j-m)T,  -  0]  -  cost  (23t-u,  )(t-mT,  )  -  0] 
lb  lb 

23t  -  u. 


j-2  cost  (23t-u,  )(i+l  )T,  -  0]  -  cost  (23t-u,  )iT,  -  0] 

,  V  j  _ 1  _ b  1  b 

L  i  23t  -  u 

i  =  j-m  1 


^j-l 


cost  (2j3t-o>1  )t  -  0]  -  cosl(23t-«1)(j-l)Tb  -  0] 
_____ 


Hence  the  output  of  the  envelope  detector  due  to  the  binary  PSK  signal  is 


given  by 


1  2,  ,  2,  .  1/2 
T(t)  =  2  [  gc  U)  +  gs  U)  . 


(m*  1  m*  1  r 

I  I  aiak  sint^j-  0)  sin(^k~  0) 

1=0  k=0  L 


+  cos(0.-  0)  cos (<p  -  0) 
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120 


t  12(m+l)  "  2  t  d.  d. 

3t-o,  ]  L  j-m+i-2  j-m+i-1 

II  1=1 


m+1  m+1 


I  I  aiVos(*i  -  V 


1/2 


1=0  k  =0 
1  *  k 


In  (28)-(30),  m+1  i  j  i  M  and  (j-I)Tb<  t  £  jTb<  Also, 


a  =  -  d  . 
o  j-m-1 


a.  s  d .  .  0  -  d .  ;  lsism 

l  j-m+i-2  j-m+i-1 


a  .  =  d .  , 
m+1  j-1 


(30) 


(31) 

(32) 

(33) 


<pQ  =  (20t  -  «  )(t  -  mTb)  (34) 

s  (20t  -  «  )(j-m-l+i)Tb  ;  1  s  i  s  m  (35) 

and 

<t>  ,  s  (20t  -  «,)t  .  (36) 

m+1  1 

Eqs.(28)-(36)  can  now  be  used  in  (15)  to  determine  the  conditional 
probability  of  detection  in  the  presence  of  a  particular  pattern  of  the 
data  sequence.  The  probability  of  detection  for  a  single  sample  can  then  be 
obtained  by  averaging  over  all  the  possible  data  sequences.  Similarly,  if 
the  carrier  frequency  is  unknown  and  hence  multiple  observations  are  used 
to  detect  the  presence  of  the  binary  PSK  data  signal,  (18)  can  be  used  to 
obtain  the  conditional  probability  of  detection  and  then  the  overall 
probability  of  correct  decision  is  given  by 
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pd  =  e  i-na-pd)  . 

1=1  i 


where  the  expectation  is  taken  with  respect  to  all  the  possible  outcomes  of 
data  bits  during  the  observation  interval.  The  maximum  number  of 
independent  samples,  N,  is  again  given  by  (25). 


III.  ANALYSIS  OF  RADIOMETER 

Consider  the  radiometer  illustrated  in  Fig.3,  where  the  input  X(t) 
may  either  be  signal  plus  noise  (H^)  or  noise  alone  (H^).  As  in  Section 
II,  the  noise  component  n^lt)  is  simply  modeled  as  additive  white  Gaussian 
noise  with  double-sided  spectral  density  t/q/2.  The  ban,  filter  is 

assumed  to  be  an  ideal  rectangular  filter  with  center  frequency  f  and 
bandwidth  2B.  The  output  of  the  bandpass  filter  r^t)  is  given  by 


r^t)  = 


y(t)  +  n(t)  ;  , 


where  y(t)  is  the  output  component  due  to  the  signal  f(t),  and  n(t)  is  the 
noise  output.  The  waveform  r ^ ( t )  is  then  used  as  an  input  to  a  square-law 
device,  and  the  output  r^ft)  is  integrated  during  the  observation  interval 
(0,  T)  to  obtain  g(T).  The  receiver  will  decide,  based  on  the  sample  value 
g(T),  whether  or  is  chosen. 


A.  Tone  Signal 

If  the  input  signal  f(t)  is  a  pure  tone,  the  performance  of  a 
radiometer  is  well-known  (l!(31.  However,  the  analysis  is  outlined  briefly 


below  because  it  leads  directly  into  the  succeeding  analysis  on  the 
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detectability  of  a  binary  PSK  signal.  Because  f(t)  is  a  tone,  i 
through  the  bandpass  filter  without  any  distortion.  Hence, 

y(t)  =  A  cos^  Ujt  +  8  j  . 

The  Gaussian  noise  n(t)  in  (38)  can  be  represented  as 


n(t)  =  n  (t)  cos 
c 


art  +  0  ]  -  n  (t)  sin  art  +  6 
1  J  s  (1 


Note  that  both  n^ft)  and  ns(t)  are  independent,  zero-mean 
processes  with  variance  VgB.  When  the  signal  is  present  and  a 
g(T)  is  approximately  given  by 

g(T)  =  |  Jj  [a  +  nc(t)j  +  ns2(t)  |  dt  . 


Since  A  +  n  (t)  and  n  (t)  are  bandlimited,  they  can  be 
c  s  J 

represented  by  their  sample  values  (21,  i.e., 


and 


A  + 


n  (t)  = 
c 


A  +  n  (i/2B) 
c 


sinc(2Bt  -  i) 


n  (t)  =  V  n  (i/2B)  sinc(2Bt  -  i)  , 

S  La  S 

1  =  -CD 

where 

.  ,  .  _  sin  (  ttx  ) 

sinc(x)  =  — -  . 

TTX 


Using  (42)-(44)  and  the  following  approximation  [  1 1, 


can  pass 

(39) 

(40) 

Gaussian 
>  T  »  271, 

(41) 

completely 

(42) 

(43) 

(44) 


sinc(2Bt  -  i)  sinc(2Bt  -  j)  dt 


I 


1 

2B 

0 


,  i  =  j  and  1  s  i  s  2BT 
,  otherwise  , 


(41)  can  be  rewritten  as  [1] 


.  2BT  /  r  -|  2  \ 

g(T)  =  4g"  [  {  [A  +  nc(i/2B)j  +  nsZ(i/2B)  |  . 


(45) 


(46) 


Note  that  n  (t)  and  n  (t+i/2B)  are  uncorrelated  for  integer  values  of  i, 
c  c 

and  thus  are  independent.  The  same  statement  is  also  true  for  ng(t). 
Hence,  the  conditional  expected  value  and  the  conditional  variance  of  g(T) 
are  given,  respectively,  by 


4-  [  a2  *  4%  ] 


(47) 


and 


V  [  A" +  2%  • 


(48) 


By  invoking  the  central  limit  theorem,  we  have 

v  -  2BT 
vr2BT 

and 


o  ~  f  v  -  2BT  -  ABT  ’I 
PH  =  <P - -  — - - 

^  V  2BT  +  2ABT  > 


where  v  is  the  threshold  of  the  test. 


X  = 


A  2  /2 

V 


is  signal-to-noise  power  ratio  and 


(49) 


(50) 


(51) 
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X 

exp(-y2/2)  dy  .  (52) 

“CO 

B.  Binary  PSK  Signal 

Suppose  that  f(t)  is  the  binary  PSK  signal  given  in  (26).  The 
impulse  response  of  the  rectangular  bandpass  filter  is  given  by 


h(t)  =  4B  sinc(2Bt)  cos(cjQt)  . 


(53) 


Hence,  the  output  y(t)  of  the  bandpass  filter  due  to  the  signal  f(t)  may  be 
obtained  by 


y(t)  = 


CD 

f(t)  h(t  -  x)  dx 

-00 


=  A  u(t)  cos(ojQt  +  0)  , 


(54) 


where 


u(t)  = 


i  m-i  r 

.-Ltd. 

n  L  l 
i  =o  L 


Si( 2tiB[ t— iT.  ] )  -  Si(27iB[t-(i+l)T.  ]) 
b  b 


(55) 


and 


Si(x)  = 


sin(y) 


dy  . 


(56) 


Under  hypothesis  H^,,  we  may  write  g(T),  analogous  to  (46),  as 


g(T) 


,  2BT  I  r 

iri  {  Au(l 

1  =  1  V  L 


Au(i/2B)  +  n  <i/2B)l  +  n  Z(i/2B)  \  . 

c  J  s  J 


(57) 


Hence,  the  conditional  expectation  and  the  conditional  variance  of  g(T) 
under  H^,  given  all  the  data  bits  during  the  time  interval  (0,  T),  can  be 
expressed  as 
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IV  NUMERICAL  RESULTS 
A.  Tone  Signal 

The  performance  of  each  receiver  will  be  illustrated  by  presenting 
its  receiver  operating  characteristic  (ROC),  which  shows  probability  of 
detection  versus  probability  of  false  alarm  [8].  Suppose  that  the 
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frequency  of  the  tone  signal  is  known  and  the  sampling  time  is  perfect, 


that  is,  tg  =  u^/23.  If  Hj  is  true,  then  the  output  of  the  envelope 
detector  due  to  fit)  will  have  its  maximum  value, 

e  (u  /23)  =  AT  /4  .  163, 


When  the  value  of  P^.  is  specified,  can  be  expressed  as 


Pd  =  q|  v2ABT^,  /  -21n(Pfa)  , 


(64) 


where  A  is  the  signal-to-noise  power  ratio  given  by  (51).  Figs.4a-4c 

illustrate  the  ROC  of  the  two  receivers  for  various  values  of  A.  The 
parameters  are  chosen  as  follows:  BT  =  600,  T/T  =  2,  A  =  -I5dB,  -20dB,  and 
-25dB.  When  the  values  of  A  increase,  the  curves  of  both  receivers  move 
toward  the  upper-left  corner,  as  expected.  However,  in  all  cases 

considered,  the  compressive  : — eiver  has  a  better  performance.  That  is, 
for  a  given  value  of  Pj.  ,  the  compressive  receiver  has  a  higher  value  of 

detection  probability.  The  result  is  quite  reasonable  because,  roughly 

speaking,  a  compressive  receiver  is  equivalent  to  a  narrowband  filter  bank 
that  has  a  higher  output  of  signal-to-noise  ratio  than  does  a  wideband 

radiometer. 


Ir  the  sampling  time  is  in  error  by 


of  main  lobe  in  (13),  P,  in  (15)  becomes 

d 


103T 


,  which  is  57.  of  the  width 


Pd  -  0 


f  J~ 2ABT1  |  sinc(  1/10)  | ,  / -21n( P f  ^ ) 


(65) 


Figs.5a-5b  show  that  the  effect  of  this  error  is  negligible. 

Suppose  now  the  frequency  .3  unknown  so  that  multiple  observations 


are  used  to  detect  the  presence  of  the  tone  signal 


The  ROC  is  shown  in 


Figs.6a-6c  for  the  same  parameters  as  those  of  Fig. 4.  Note  that  two  curves 

are  displayed  for  the  multiple  sample  case  in  each  figure;  they  correspond 

to  the  best  case  (where  one  of  the  samples  is  located  at  t  =  and  the 

worst  case  (where  one  of  the  samples  is  located  at  t  =  c^/2f3  +  A/2).  Note 

also  that  the  samples  are  uniformly  spaced  at  an  interval  of 

A  =  jt/0T  .  (66) 

Even  under  the  best  situation,  the  curves  in  Figs.6a-6c  show  that  the 

compressive  receiver  performs  worse  for  multiple  observations.  The  reason 

for  this  is  that,  if  the  overall  probability  of  false  alarm,  P_  ,  is  fixed, 

r  A 

the  value  of  P^.  for  an  individual  sample  must  decrease.  Hence,  for  a 
fixed  value  of  A,  the  threshold,  y,  of  the  test  must  be  higher  for  multiple 

observations,  and  this  decreases  the  probability  of  detection,  P  ,  of  each 

sample. 

B.  Binary  PSK  signal 

If  we  assume  that  the  carrier  frequency  is  known  and  the  sampling 

time  of  the  compressive  receiver  is  t  =  ca,/2/3,  g  (t  )  and  g  (t  )  in 

s  1  &c  s  s  s 

(28M29)  can  be  simplified  as 

*clts>  -  5  cos(6){d  ^(jT  -  V  .f  d  Tb  *  d  (t  -  u-iny} 

V  J  1  =  J -m  J  > 

AT1  - 

=  —  I  cos(0)  (67) 

and 

AT1  - 

g  (t  )  =  1  sin(9)  ,  (68) 

S  S  Z 

where 
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d(x)  dr 


(69) 


t 

a 

t  -T 

a  1 

is  the  time  average  of  the  data  bits  during  the  time  period  T^.  Note  that 
T  is  a  random  variable,  since  it  depends  upon  the  data  bits  in  the  window. 
Substituting  (67)-(68)  into  (30),  we  have 

e.(t  )  =  AT,  III /4  .  (70) 

f  s  1 1  1 

When  the  value  of  P^a  is  specified, 

Pd(|I|)  =  q(  /  2ABT1  |I j ,  /  -21n(Pfa)  1  (71) 


is  the  conditional  probability  of  detection,  given  the  value  of  |T|.  The 
unconditional  probability  of  detection  for  the  compressive  receiver  can 
then  be  obtained  by  taking  the  expectation  with  respect  to  III,  i.e., 


pd'E(pd(l'l))  • 


(72) 


Although  (71 )— (72)  are  derived  for  mT^  =  T^,  where  m  is  an  integer,  they 
are,  in  fact,  also  valid  if  m  is  a  real  number.  Assuming  k  is  an  integer 
and  k  i  m  <  k+1,  we  may  express  (72)  as 


P  =  (1  -  ^)  P  (1)  +  ^ 
d  2  d  2 


Pd(|l  -  2x|)  dx  ;  k  =  0 


(73) 


and 


pd  =  [“=«  (  i )  [(m  - kl  pd(l‘  ‘  Zi/ml’ 


m-k 


Pd(  1 1  -  2i/m  -  2x/m|)dx 


i  i « ( v )  [<* 


■»  1  -  m)  P  1 1  -  2i/m  | ) 
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m-k 


P^(  1 1  -  2i/m  -  2x/m  j  )dx 


]' 


k  a  1  , 


(74) 


where 

( j )  *  jT-reWn  •  1751 

and  the  conditional  detection  probability  P  ( • )  is  again  given  by  (71). 
The  corresponding  probabilities  for  the  radiometer  have  been  derived  in 
Section  III.  Hence,  we  can  compare  their  performance  by  plotting  the  ROC 
for  both  receivers.  The  parameters  in  Figs.7a-7c  are  chosen  as  follows:  BT 
=  600,  T/T  =  2,  T^T^  =  4,  X  =  -lOdB,  -15dB  and  -20dB.  Unlike  the  previous 
case,  the  performance  of  the  compressive  receiver  is  not  always  superior 
to  that  of  the  radiometer.  This  is  especially  true  when  P^a  is  large  and  X 
is  high.  However,  in  practice,  we  would  typically  choose  a  receiver  which 
has  a  better  performance  for  small  values  of  X  and  P^a  (say,  P^  = 
1.0E-06).  Under  this  situation,  the  compressive  receiver  outperforms  the 

radiometer. 

If  we  compare  the  curves  of  the  compressive  receiver  in  Fig. 7  with 

the  corresponding  curves  in  Fig.  4,  we  see  that  the  performance  for  the 

binary  PSK  data  signal  is  much  worse  than  it  is  for  the  tone  signal.  The 

degradation  is  obviously  from  the  term  I.  Since  each  data  bit  d.  is  equal 

to  +1  or  -1  with  probability  0.5,  the  absolute  value  of  the  time  averap'' 

these  data  bits  is  usually  smaller  than  1.  On  the  other  hand,  the  KG^  of 

the  radiometer  is  almost  the  same  for  both  signals.  In  iact,  if  BT  »  1, 

the  intersymbol  interference  components  of  u(t)  in  (15)  are  negligible  and 
2 

u  (t)  is  approximately  equal  to  1.  Hence,  P^  in  (61)  is  insensitive  to 

values  of  data  sequence  (d.)  and  (62)  reduces  to  (50). 
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In  Figs.8a-8c,  the  ROC  for  both  receivers  are  displayed.  The 


parameters  are  chosen  to  be  the  same  as  those  of  Fig. 7  except  that  M  =  8. 
The  curves  in  Figs. 7  and  8  are  very  similar.  We  see  that  the  detection 
probability  of  the  radiometer  is,  as  expected,  almost  invariant  to  the 
change  of  M,  but  the  compressive  receiver  has  a  worse  performance. 

Suppose  now  the  carrier  frequency  of  the  binary  PSK  signal  is 
unknown  and  N  samples  are  used  to  detect  the  presence  of  the  signal. 
Similar  to  the  case  of  the  tone  signal,  the  performance  of  the  compressive 
receiver  depends  upon  the  values  of  the  sampling  time.  Note  that  one  of 
the  N  samples  will  be  obtained  from  the  interval  [u^/2/3,  Wj/2/3  +  A),  where 
A  is  again  given  by  (66).  The  effect  of  the  offset  from  w^/2/S  is 
illustrated  in  Fig. 9,  for  the  case  that  the  carrier  frequency  is  at  the 
center  of  the  band.  The  values  of  the  legend  represent  the  offset  in  terms 
of  A.  It  is  clear  that  the  curve  with  zero  offset  represents  the  best  case 
and  the  one  with  0.5A  offset  corresponds  to  the  worst  case. 

The  results  of  the  compressive  receiver  with  multiple  samples  (N  = 
600)  are  also  shown  in  Figs.7a-7c.  For  A  =  -lOdB,  the  compressive  receiver 
with  multiple  samples  performs  better  the  other  two  receivers.  When  A  = 
-15dB,  the  compressive  receiver  has  a  performance  somewhere  between  that  of 
the  compressive  receiver  with  a  signal  sample  and  the  radiometer.  If  A  = 
-20d&\  the  compressive  receiver  with  multiple  samples  has  the  worst 
performance.  In  other  words  the  compressive  receiver  with  multiple  samples 
appears  to  perform  better  than  one  using  just  a  single  sample  for  high 
signal-to-noise  power  ratios.  This  is  interesting,  because  the  frequency 
of  the  signal  is  assumed  known  to  the  receiver  in  the  single  sample  case, 
and  is  assumed  unknown  in  the  multiple  sample  case.  The  reason  for  this  is 
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that  the  samples  in  the  sidelobes,  when  A  is  high,  are  useful  in  detecting 
the  presence  of  the  signal. 


V.  CONCLUSIONS 

The  performance  analysis  of  a  compressive  receiver  when  used  to 
detect  the  presence  of  either  a  pure  tone  or  a  binary  PSK  signal  has  been 
presented.  The  receiver  operating  characteristic  of  the  compressive 
receiver  was  compared  with  that  of  a  wideband  radiometer.  If  the  carrier 
frequency  of  the  signal  is  known,  the  compressive  receiver  uses  a  single 
sample  to  detect  the  presence  of  the  signal.  Otherwise,  multiple  samples 
are  used  to  determine  whether  or  is  true. 

When  a  tone  signal  is  to  be  detected  and  the  frequency  is  known,  the 
compressive  receiver  performs  better  than  the  radiometer.  We  have 
demonstrated  that  the  affect  of  a  small  offset  of  the  sampling  time  of  the 
compressive  receiver  is  negligible.  Also,  it  appears  that  multiple 

observations  cannot  improve  the  detectability  of  the  compressive  receiver 
for  a  tone  signal. 

For  a  binary  PSK  signal,  even  if  the  frequency  is  known  to  the 

compressive  receiver,  the  performance  of  the  compressive  receiver  is  not 

always  superior  to  that  of  the  radiometer.  However,  under  the  practical 
conditions  of  small  values  of  signal-to-noise  power  ratio  and  false-alarm 
probability,  the  compressive  receiver  outperforms  the  radiometer.  For  the 
case  of  multiple  samples,  the  compressive  receiver  has  better  performance 

than  the  radiometer  when  A  is  high,  and  also  outperforms  the  single-sample 
compressive  receiver.  For  small  values  of  A,  the  single-sample  compressive 
receiver,  operating  with  the  optimum  sampling  time,  is  superior  to  one 
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using  multiple  samples.  However,  if  the  carrier  frequency  is  unknown,  this 
latter  case  is  of  no  practical  interest. 
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Fig. 8b  ROC  for  BPSK  Signal 
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Fig. 8b  ROC  for  BPSK  Signal 
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Ac  ?  t  r  3  :  t 

A  method  ftr  applying  spread  ;pe::rrr, 
•.e:r.r.:rjes  ::  optical  ccmmunicat ior.  is 
creser.ted.  Tr.e  interference  suppression 
capability  of  spread  spectrum  systems  os 

domain  processing  and  a  design  for  spectral 
coding  an  optical  frequencies  is  giver..  The 
encoding  system  is  cased  on  optical  pulse 
compression  and  Shaping.  Several  possible 
implementations  of  this  system  are 
suggested,  and  applications  to  fiber  optics, 
laser  radar,  free  space  optical 
communications ,  and  cor.er  systems  are 
diseased . 


Introduction  -  Spread  Spectrum.  Systems  and 
Trar. sform.  Domain  Processing 

Spread  spectrum,  communication  consists 
of  transmitting  a  given  signal  by  modulating 
tr.e  information  with  a  large  bandwidth, 
coded  waveform  such  as  a  PN  sequence  [1-4; . 
Tr.e  transmitted  signal  occupies  a  bandwidth, 
much  larger  than  tr.e  information  bandwidth. 
Suer,  systems  posses  a  number  of  special 
properties  w.hicn  distinguish  them,  from 
narrtwcar.d  communication  techniques.  A 
primary  advantage  cf  such,  systems  is 
resistance  to  ;  anting  and  interference.  A 
bread  spectral  bandwidth,  signal  13  more 

ficult  to  distinguish  from  ambient  noise, 
wr.ic.h  adds  tc  t.ne  security  of  the  channel. 

Spread  spectrum  techniques  have  net 
peer,  utilized  m  optical  communication 
systems,  despite  their  increasing  popularity 
and  inherently  wide  bar.dwidtr.s,  due  to  a 
la:*  cf  effective  modulation  and  coding 
methods  available  at  optical  frequencies. 

We  propose  several  methods  for  incorporating 
spread  spectrum  techniques  into  optical 

A  typical  direct  sequence  spread 
spectrum,  communications  system,  is  sr.owr.  m 
Fig.  1.  hr.  .nfcrmaticr.  sequence,  S(ti,  is 
tcou.aced  cy  a  :;oe  sequence,  r(t)  .  7-e 


•This  wer*  r.a3  peer,  supe  .r-or  m  cart  :v 


communications  channel  by  interference, 

Tit),  and  additive,  almost  white  Gaussian 
noise,  n(t) .  The  corrupted  signal  is 
recovered  by  a  matched  filter  containing  the 
code  sequence . 

The  ability  cf  a  spread  spectrum 
system,  to  resist  jamming  is  determined  by 
the  processing  gain,  which  in  turn  is  giver, 
cy  the  ratio  of  transmission  bandwidth  to 
data  bandwidth.  Large  processing  gams 
provide  a  high  degree  of  jamming  immunity. 
Since  processing  gain  cannot  be  increased 
rndef initely,  it  is  desirable  to  supplement 
the  jamming  resistance.  This  has  led  to  the 
use  of  transform  domain  processing 
techniques  [4-10], 

A  transform  domain  receiver  is  shown 
in  Fig.  2.  The  received  signal  is  s(t)  - 
S(t)c(t),  plus  channel  interference  and 
noise.  Filtering  by  the  transfer  function 
H(Ci)i  is  performed  by  multiplication  followed 
by  inverse  transformation.  This  real  time 
frequency  domain  multiplication  ha3  beer, 
demonstrated  both  theoretically  and 
experimentally  (1,2].  An  alternate  receiver 
implementation  replaces  the  matched  filter 
by  multiplication  with  the  signal  spectra  in 
the  transform  domain.  Talcing  the  inverse 
Fourier  transform  produces  a  result 
equivalent  to  Fig.  2. 

Transform  domain  processing  techniques 
effectively  suppress  narrow  band  jammers  in 
a  spread  spectrum  system.  The  jammer  may  be 
removed  by  the  system  illustrated  in  Fig.  3. 
Input  consists  of  the  code  and  jammer  on  an 
RF  carrier;  a  high  power,  narrow  band  jammer 
appears  as  an  impulse  in  the  transform 
domain.  A  gating  function  removes  the 
portion  of  the  spectrum  containing  the 
;ammer.  The  gate  output  is  the  PN  code; 
since  the  code  has  a  large  bandwidth,  the 
notch  filter  has  not  seriously  degraded  the 
signal  spectrum.  Correlation  is  performed 
by  multiplication  in  the  transform  domain, 
followed  by  ar.  inverse  Fourier  transform. 


Ipcical  signal  processing  techniques 
ire  juiceo  tc  applications  m  transform 
terrain  processing.  We  shall  consider 
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where  c  i3  the  speed  of  light, 
the  optical  pulse  becomes 


he  phase  of 


Chirp  Transform,  which  car.  be  used  to 
implement  a  real-time  Fourier  transform.. 

The  Chirp  transform  is  illustrated  by 
Fig.  4(a).  An  arbitrary  signal  f(t)  is 
multiplied  by  a  down-chirp,  then  passes 
through  a  linear  system  whose  impulse 
response  is  an  up-chirp.  The  result  is 
multiplied  by  down-chirp;  the  output  is  the 
Fourier  transform  of  the  input.  This  is 
aisc  known  as  the  mu  It  iply-cor.vol  ve-mult  iply 
(v;CM)  algorithm;  the  system  of  Fig.  -5  (b) 
performs  the  same  operation  by  a  dual 
process  called  the  convolve-muit iply- 
ecr.vc.ve  (CMC)  algorithm. 

Either  algorithm  can  be  used  in  optics 
to  cerform  spatial  Fourier  transforms  of 
ccticai  signals.  The  MCM  algorithm  may  also 
be  applied  to  optical  pulse  com.press.on 
s vs terns  This  process  can  be  realized  in 

the  time  domain,  which  forms  the  basi3  for 
me  use  of  spread  spectrum  systems  at 
optical  frequencies. 


The  basic  configuration  fcr  optical 
pulse  compression  and  shaping  is  shown  in 
Fig.  :,  after  the  treatment  of  Ref.  (11,12;. 
N'ote  that  Fig.  5  is  an  optical 
imp lemer.ta tion  cf  Fig.  4(a),  w.nere  f(t)  is 
me  pulse  to  be  compressed.  This  pulse 
com.cressicn  is  achieved  by  inducing  a  chirp, 
or  linear  frequency  sweep,  or  an  optical 
culse  and  subsequently  rephasing  the  chirped 
frequency  components.  Although  there  are 
several  means  cf  obtaining  chirped  optical 
pulses,  a  single  mode  optical  fiber  induces 
uniform  frequency  modulation  across  the 
entire  pulse  profile  The  frequency  chirp 
is  generated  by  self-phase  modulation,  which 
arises  from  the  interaction  of  the 
propagating  light  and  the  intensity 
dependent  portion  of  the  fiber's  refractive 
index  r  t  3 ,  .  It  is  then  necessary  to  re- 
cnase  the  spectral  components  to  compress 
me  pulse  in  time.  The  system  in.  Fig.  5 
uses  a  diffraction,  grating  pair  as  a 
dispersive  delay  line  (14].  This  system  can 
to  or.:: ie  an  optical  pulse  in 
rne  frequency  den.  a  m  . 

A  chirped  optical  pulse  can  be 
produced  by  fhe  nonlinear  process  c-  sel*~ 
pr. as-  .-.m-ulation  <S?M)  .  Tiie  fiber's 
ref:  mve  index  is  given  by 

n  =  no  — 


wr.ere  ' 


is  the  intensity  profile  of  the 
r.T  is  a  positive  materia. 
y  .  The  propagation  constant  is 


(2) 


e  =  wot  -  -  AI{t)  (3) 

c 

where  A  represents  a  collection  of  terms  and 
z  is  the  propagation  distance.  The 
instantaneous  frequency  is  thus  proportional 
to  the  negative  time  derivative  of  the 
intensity  profile, 

-  -  di  =  **  -  4  w  w 


and  the  properties  of  the  resulting  chirp 
depend  on  the  time-varying  intensity. 

However,  the  chirp  produced  by  I ?.M. 
alone  is  not  linear  over  the  full  intensity 
profile.  The  linearity  of  the  chirp  can  be 
improved  by  the  effect  of  positive  group 
velocity  dispersion  (GVD)  in  the  fiber. 

This  effect  is  calculated  by  expanding  the 
propagation  constant,  k(w),  about  the  center 
frequency  Uq  ;  the  relevant  term  is 


0  = 


k2L 


(5) 


The  combined  effect  of  SPM  and  GVD  13  called 
dispersive  self-phase  modulation  (DSPM) ;  it 
produces  an  approximate  square,  linearly  up- 
chirped  pulse  from  a  single  frequency  input 
of  nonuniform  intensify. 

As  mentioned  before,  the  system  of 
Fig.  5  implements  the  following  relation: 


(/(0«+,v,t*)  *  (6) 

where  f(t)  is  the  modified  light  pulse,  F  (0)1 
is  its  Fourier  transform,  U)  •  2hf  is  the 
angular  frequency,  and  (3  is  a  chirping 
factor  to  be  determined. 

Convolution  is  equivalent  to  passing 
the  chirped  signal  through  a  linear  system 
with  a  chirp  as  its  impulse  response.  In 
the  frequency  domain,  this  is  the  same  as 
multiplying  by  the  Fourier  transform  cf  the 
chirp  signal. 


H(f)  = 


V  0 


(') 


The  grating  pair  acts  as  a  delay  line; 
its  frequency  response  is 

Hl(f)=e^T  (S) 


where  T  is  the  delay  for  a  given  frequency 
component.  By  equating  the  exponents  cf  eq 
(7)  and  (9),  the  grating  pair  with  a  chirp 
input  will  complete  the  realization  cf  the 
left  .-.ana  side  cf  eq.  ii)  creviced  t.-.-i 
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Ar  =  t  A//'/?  (9) 

where  At  is  the  difference  in  delay  time 
between  two  spectral  components  with 
frequency  difference  Af.  This  is  a  standard 
relation  for  diffraction  gratings  (14 ], 
which  indicates  that  the  grating  transit 
time  is  inversely  proportional  to  frequency. 
It  is  this  effect  which  makes  temporal  pulse 
compression  with  a  grating  pair  possible. 

The  chirp  parameter  for  a  grating  pair 
is  giver,  by  (14( 

3  =  tt(A /)-  d/bT  (10) 

wr.ere  c  is  t r.e  grating  constant  *  1/ nurr.be r 
: :  lines  per  rrrr.,  n  is  the  grating 
separation#  and  T  is  a  const an c  determined 
ry  the  angle  of  incidence.  Thus,  the  system 
i  Fig.  5  realizes  eq.  (6);  the  analogous 
functions  for  each  component  are  shown  in 
Fig.  r.  After  one  pass  through  the  grating 
pair  the  optical  signal,  g^(t),  is  given  by 

9i(t)  =  F{3t,  r)  e-^'5  (11) 

The  first  pass  through  the  gratings  has 
taxer.  tr.e  Fourier  transform  of  the  original 
ut  signal,  because  the  grating  pair  input 
the  opposite  chirp  as  the  impulse 
response,  hit),  of  the  grating  pair.  The 
signal  of  eq.(6)  .now  possesses  a  chirp  cf 
tee  same  sign  as  hit);  if  this  signal 
reflects  from,  the  mirror  and  traverses  the 
grating  pair  a  second  time,  no  further 
Fourier  transformation  can  occur.  This 
second  pass  inrough  the  gratings  doubles  the 
or.irp  factor;  tr.e  final  output  of  the 
system,  g; ( t ) ,  is 

g2(t)  =  F(3t;z)  e~l23t'  (12) 

The  system  has  performed  pulse 
compression  as  an  y.Cy.  system.  After  one 
pass  through  the  gratings,  the  Fourier 
spectra  of  the  signal  is  spatially  separated 
m  a  plane.  3y  placing  seme  form  of 
transmission  mas*  in  this  plane,  it  is 
possible  to  modulate  the  frequency 
components  of  the  optical  signal.  If  a 
t r ar.sm.iss ion  masx  function  Mlf)  is  placed  in 
this  plane  as  indicated  by  Fig.  7,  then  the 
final  output  will  be 

9;(t)  =  F(3i  m)  A/  (3t;ir)  e~l2jr  (13) 

heritage  et.al.  lo(  nave  experimented  with 
simp.e  cwo-dimensicr.al  amplitude  and  phase 
masxs,  as  well  as  spectral  windows  to 
eliminate  nonlinear  ber.avior  at  the  euges  of 
tr.e  pulse  (Is)  .  We  propose  to  use  this 
principle  to  achieve  spread  spectrum  coding 
of  tr.e  optical  signals. 


Note  that  in  place  of  a  grating  pair, 
a  fiber  optic  delay  line  with  negative  3VD 
could  be  used  if  only  pulse  compression  was 
desired.  However,  the  grating  pair  is  more 
suitable  for  our  purposes,  as  it  provides  a 
means  to  modify  the  spectral  components . 

Spread  Scectrum  Techr.icues  ir.  On-  •  -al 
Co.-rmunicat  ion 

There  are  many  possible  designs  for  an 
optical  spread  spectrum  communication  system 
based  on  the  pulse  compression  architecture 
described  previously.  Fig.  7  illustrates 
how  the  chirped  light  may  be  modulated  with 
binary  data;  the  appropriate  pn  code  must 
now  be  impressed  upon  this  data.  Cne 
implementation  involves  placing  a 
transmission-type  mas*  between  the  grating 
pair  and  the  feedback  mirror,  as  shown  m 
Fig.  7.  Since  the  frequency  components  are 
spatially  distributed  in  this  plane,  it  is 
possible  to  perform  both  amplitude  and 
phase  modulation  on  the  optical  signal.  The 
optical  signal  passes  through  the 
transmission  mas*  twice;  the  second  pass  oar. 
only  be  neglected  if  a  binary  amplitude  m.asx 
is  used  (consisting  of  either  opaque  or 
transparent  pixels) .  If  more  complicated 
masks  are  required,  as  in  the  case  of  spread 
spectrum,  systems,  then  the  square  root  of 
the  desired  amplitude  function  must  be 
implemented  on  the  mask.  Any  type  of  phase 
coding  must  account  for  the  double  phase 
delay  incurred  by  a  second  pass  through  the 
mas*.  The  optical  signal  at  the  encoding 
plane  is  the  Fourier  transform  of  the  input 
pulse;  thus,  the  encoding  mas*  must  be  the 
Fourier  transform  of  the  desired  code.  The 
output  of  the  system  is  then  giver,  by 
eq.  (8)  . 

The  encoded  optical  pulse  is 
transmitted  along  a  fiber  optic  link.  v;e 
assume  that  any  further  chirping  due  to  this 
transmission  may  be  accounted  for.  An 
optical  receiver  and  decoding  scheme  for 
this  signal  is  shown  in  Fig.  3.  The  optical 
signal  is  passed  through  a  grating  pair, 
which  spatially  separates  the  frequency 
components  without  affecting  the  encoded 
signal.  The  grating  pair  is  now  separated 
by  twice  the  distance  provided  at  the 
encoder,  to  account  for  the  doubled  chirp 
rate.  This  spatially  dispersed  signal  is 
then  passed  through  the  complex  conjugate  cf 
the  encoding  mask,  M*(f)  .  Multiplication  ir. 
the  transform  domain  is  equivalent  to 
correlation  in  the  time  domain.  If  the 
decoding  mas*  matches  the  signal  modulation, 
a  correlation  pea*  will  be  observed. 
Otherwise,  the  output  will  resemble  random 
noise,  since  the  crcss-correlation  of  two 
different  codes  is  near  rerc.  The  optica, 
signal  passes  through  another  grating  pair 
to  compress  the  spectra  befere  oetecticn. 


The  correlation  receiver  design  is  greatly 
simplified  ir.  tr.e  transform  domain.  Spread 
spectrum  ceding  techniques  at  optical 
frequencies  can  be  implemented  in  this  way. 

An  alternative  approach  eliminates  the 
second  pass  through  the  grating  pair  at  both 
the  transmitter  and  receiver.  This  second 
pass  serves  to  recompress  the  optical 
spectra  spatially  prior  to  transmission,  and 
also  doubles  the  chirp  rate.  In  the  system 
cf  rig.  9,  a  ler.s  is  placed  behind  the 
transmission  mask  to  spatially  recombine  the 
spectra.  The  transmitted  signal  retains  its 
original  chirp  rate,  so  that  the  grating 
spacing  at  both  the  transmitter  and  receiver 
may  be  the  same.  The  receiver  design  is 
similar  to  rig.  9;  a  grating  pair  spatially 
disperses  tr.e  optical  pulse,  which  passes 
through  a  decoding  mas*.  A  lens  replaces 
the  second  grating  pair  to  compress  the 
spatially  distributed  pulse  for  optical 
detection.  This  system  would  require  a 
nearly  "ideal"  lens  at  both  the  transmitter 
and  receiver  -  a  lens  possessing  a  large 
field  of  view  and  practically  no  chromatic 
dispersion.  Compound  lenses  of  this  type 
are  availacle,  although  they  are  much  more 
expensive  than  high  quality  diffraction 
gratings.  It  has  been  noted  that  a  single 
pass  through  the  grating  pair  produces  a 
slightly  elliptical  beam  profile  [14];  this 
is  corrected  by  a  second  pass  through  the 
grating  pair. 

These  systems  transmit  coded  optical 
pulses  in  the  transform  domain,  and  decode 
them  in  the  same  manner.  For  some 
applications,  we  may  transmit  the  temporal 
form  of  the  ceded  pulses.  This  requires  an 
inverse  Fourier  transform  at  the 
transmitter.  The  mirror  behind  the  mask  has 
been  replaced  by  a  phase  conjugate  mirror! 
I'pcn  reflection  from,  such  a  mirror,  the 
phase  of  the  encoded  pulse  will  be  reversed; 
after  a  second  pass  through  the  gratings, 
the  modulated  pulse  is  returned  to  the  time 
domain.  The  received  pulse  is  passed 
through  a  grating  pair,  which  performs  a 
Fourier  transform,  and  spatially  distributes 
the  frequency  components.  A  decoding  mask 
is  used,  and  the  decoded  signal  is  spatially 
compressed  for  detection  by  either  another 
grating  pair  or  a  lens. 

Another  system  for  transmission  in  the 
time  domain  is  shewn  in  Fig.  10.  The 
encoded  pulse  is  spatially  recombined  by  a 
lens  and  focussed  into  another  length  of 
fiber  optic  cable.  The  cable  induces  an  up- 
chirp  on  the  propagating  signal.  If  a  cable 
twite  as  long  as  the  system  input  cable  is 
used,  the  dowr.-chirp  of  the  coded  signal 
will  ce  transformed  into  an  up-chirp.  This 
signal  passes  through  another  grating  pair, 
wnic.n  performs  the  inverse  Fourier 
transform.  One  pass  through  the  second 


grating  pair  would  spatially  separate  the 
spectra;  two  passes  are  required  to  obtain  a 
compacted  output.  The  received  signal  may 
be  decoded  by  the  time  domain  decoders 
presented  earlier. 

There  exists  another  method  for 
generating  optical  Dulses  of  arbitrary  shape 
which  i3  fundamentally  different  from  the 
preceding  systems  [IT],  illustrated  in  Fig. 
11.  The  pulse  shape  is  controlled  by  a 
programmable  electro-optic  modulator  at  the 
system  input.  The  optical  pulse  is 
modulated  and  amplified  prior  to  passing 
through  the  spectral  rephasing  process;  a 
translating  mirror  is  used  to  control  the 
phase  of  the  signal.  Modulating  the  puise 
before  it  enters  the  pulse  compression  stage 
controls  the  resulting  pulse  shape  and 
spectral  content  [17],  The  analysis  of  this 
system  is  similar  to  our  earlier  discussion; 
this  method  represents  an  alternative  way  to 
achieve  coding  of  the  optical  pulses.  The 
electro-optic  modulator  must  be  clocked  at  a 
rate  equal  to  the  fiber  optic  cable's  data 
rate  to  achieve  practical  modulation.  3in.ee 
the  optical  data  rate  may  approach  several 
gigabits/second,  the  processing  speed  would 
be  limited  by  the  electro-optic  modulator  ia 
few  hundred  megabits /second) .  This  is  a 
disadvantage  compared  to  systems  which 
achieve  real-time  encoding  using  transform, 
domain  techniques. 

Abdications 

There  are  many  applications  fer 
transform  domain  processing  in  tr.e  cpticai 
regime.  One  possibility  is  a  fiber  optic 
communication  system  using  spread  spectrum 
techniques.  Fiber  optics  represents  a 
secure  method  of  communication  because  cf 
its  resistance  to  electromagnetic 
interference.  The  high  band  width  cf  an. 
optical  fiber  system  may  be  fully  exploited 
to  provide  maximum  security  by  employing 
optical  encoding  techniques.  The  concept  cf 
a  free  space  laser  communication  system  has 
been  proposed;  laser  signals  suffer  from 
degradation  caused  by  pulse  spreading  or 
scattering  in  the  atmosphere.  In  order  to 
transmit  information  reliably  under  these 
oonditions,  some  form  of  coding  is  highly 
desirable . 

Another  application  is  laser  radar 
systems.  A  fundamental  problem  in 
conventional  radar  systems  is  the  decrease 
m  returning  signal  amplitude  with  range. 
Although  the  problem  can  be  overcome  cy  tr.e 
use  cf  larger  antennas  and  higher 
transmission  power  levels,  there  are 
practical  limits  to  this  approach  [13]. 

Puise  compression  techniques  combined  with 
spectral  weighting  have  been  used  ir. 
microwave  radar  systems  to  overcome  this 
problem.  By  transmitting  a  chirped  pulse. 
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it  is  possible  po  use  compression  techniques 
to  concentrate  the  energy  of  the  returning 
echo  into  a  detectable  signal.  Similar 
tec.cniques  would  be  effective  in  laser  radar 
systems,  provided  the  optical  pulse  spectra 
is  encoded. 

The  coding  cf  optical  pulses  also  has 
applications  in  fiber  oppic  code  division 
multiple  access  systems  [18].  The  scheme 
uses  the  excess  bandwidth  of  fiber  optics  to 
map  low  information  rate  electrical  or 
optical  signals  into  high  data  rate  optical 
pulse  sequences.  Single-mode  fibers  are 
well  suited  to  this  application.  The 
sequences  can  be  used  to  achieve  random, 
asynchronous  access  to  a  communication 
system,  free  of  network  control  among  many 
users.  Each  user  of  a  multipie-access 
system  is  provided  with  a  unique  code,  whose 
cross-correlation  with  other  system  codes  is 
nearly  zero.  A  set  of  Optical  Orthogonal 
Codes  has  been  developed  for  this  purpose 
[15]  .  A  coding  system,  which  operates  at 
tptical  frequencies  is  essential  to  such  a 
system . 

Because  optical  transform  domain 
processing  is  suited  to  optical  transmission 
over  long  distances,  it  possesses  several 
applications  to  outer  space  systems.  The 
Strategic  Defense  Initiative  Organization 
has  expressed  a  need  for  the  transmission  of 
high  intensity,  short  duration  laser 
pulses..  This  problem  is  analogous  to  a 
long-range  radar  system;  a  similar  solution 
involves  generating  lower  intensity  laser 
pulses  encoded  with  a  frequency  chirp.  Such 
uises  could  be  directed  to  a  target  by 
“flection  from  a  pulse  compression  system, 
producing  laser  pulses  of  shorter  duration 
and  higher  intensity.  Optical  transform 
domain  encoding  would  be  necessary  to  reduce 
the  sidelcbes  of  the  compressed  pulse. 
Smaller  scale  pulse  compression  systems  have 
provided  compression  factors  of  80  or  more 
[12] .  Systems  of  this  type  would  be 
effective  if  ground  based  lasers  were  used 
to  supply  orbital  satellites  equipped  with 
large  mirrors  and  diffraction  gracing 
compressors.  The  problems  of  high  energy 
beam  propagation  through  the  atmosphere 
could  be  minimized  by  the  use  of  optical 
transform  domain  coding. 

Another  space  application  involves  a 
proposal  by  h'ASA  to  develop  a  laser  powered 
single  stage-to-orbit  space  vehicle  by  the 
twenty-first  century  [19,20].  This  vehicle, 
cr.e  Apollo  Lightcraft,  would  receive  its 
power  from  a  remote  laser  source  and 
requires  gigawatt  optical  pulses  of  a  few 
hundred  microseconds  duration.  To  realize 
efficient  transmission  of  the  optical  power 
team,  transform  domain  encoding  could  be 


used.  Fig.  12  shows  an  artist's  rendition 
of  the  ground  based  laser. 

We  may  implement  the  optical  coding 
system  in  a  compact  format  for 
telecommunications  or  optical  computing 
applications.  The  entire  system  could  be 
realized  in  an  integrated  optics 
configuration  as  illustrated  in  Fig.  13.  A 
surface  waveguide  could  generate  the  same 
type  of  frequency  chirps  as  an  optical 
fiber;  the  same  amount  of  chirping  migr.r  be 
achieved  by  a  short  segment  of  properly 
doped  surface  waveguide.  The  waveguide  car. 
also  be  designed  as  the  delay  line  media 
between  the  two  diffraction  gratings 
Gratings  could  be  realized  using  acousto¬ 
optic  or  electro-optic  methods.  Both 
effects  can  generate  programmable 
diffraction  gratings.  In  some  materials, 
the  gratings  can  be  stored  in  the  crystal 
after  the  applied  voltage  is  removed  (toe 
photoref ractive  effect).  Both  effects  can 
exist  simultaneously  in  a  piezoelectric 
crystal  resulting  in  an  acousto-elect ro- 
optic  diffraction  grating,  which  has 
advantages  over  applying  either  effect 
separately  [21].  Since  ail  of  these  methods 
modulate  the  diffracted  light,  they  could  be 
used  to  generate  masks  for  transform  domain 
encoding.  Electro-optic  modulators  have 
been  used  [17]  to  modulate  optical  signals 
before  performing  pulse  compression.  As  tr.e 
technology  of  optical  pulse  shaping 
continues  to  develop,  the  use  of  integrated 
optics  will  become  increasingly  important. 


Conclusions 

We  have  described  a  system  for 
implementing  spread  spectrum  communication 
techniques  at  optical  frequencies.  The 
frequency  encoding  is  achieved  by  a  modified 
optical  pulse  compression  system  using 
transform  domain  processing  techniques. 

This  system  is  able  to  exploit  the  high 
bandwidth  of  optical  communicat ion  systems; 
different  optical  implementations  may  be 
realized.  Applications  include  improved 
fiber  optic  code  division  multiple  access, 
free  space  optical  communication  systems, 
-aser  raoar,  integrated  optical  systems,  an-, 
space  applications  such  as  NASA's  Apcilo 
Lightcraft  project . 
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Fig.  12  -  Artist's  conception  of  ground 
based  laser  for  space  craft 
propulsion 


Fig. 


13 


Pncectuai  design  of  an  integrate 


